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ABSTRACT 


The  aetuail  flov  conditions  vithln  the  passages  of  a  l4-inch  diaiiieter, 
four>bladed  tnq>eller  of  a  mixed  flow  compressor  have  been  investigated  and 
eaoqpared  with  design  values,  which  were  based  on  axial  symmetry  and  income 
press ible,  nonviscous  flov.  The  blades  were  designed  as  bound- vortex 
sheets,  having  a  zero  tangential  component  of  vorticity,  l.e,,  the  velocity 
distribution  in  meridional  planes  is  not  disturbed  by  the  blades.  Flow 
surveys  at  the  impeller  outlet  from  hub  to  shroud  and  from  blade  to  blade 
were  made  with  a  three-dimensional,  spherical  pitot  probe  mounted  on  the 
impeller.  This  probe  and  also  merctiry  seals  for  transferring  the  pressure 
signals  from  the  rotating  shaft  to  stationary  manometers  were  especially 
developed  for  this  investigation.  Test  runs  conducted  in  the  range  from 
1500  to  3000  rp®  8®d  from  8.7  to  44.1  cubic  feet  per  second  shoved  that 
the  shape  of  the  meridional  velocity  distribution  between  hub  and  shroud 
agrees  essentially  with  the  theory.  However,  in  the  blade-to-blade  direction 
the  distribution  curves  for  the  meridional  velocity  are  displaced  as 
much  SIS  20^  above  and  10^  below  the  average.  Vauriations  in  the  cOne 
angle  of  the  relative  velocity  are  only  in  the  order  of  2°  and  3°»  i.e.> 
the  flow  does  not  "twist"  away  much  from  surfaces  Of  revolution*  The 
median  impeller  slip  factor  was  0,57  within  the  test  range. 


The  blade-to-blade  distribution  of  the  flow  on  the  mean  stream  surface 
was  also  calculated  by  the  relaxation  method.  It  waa  found  that  the  actual 
flow  was  more  uniform  than  calc\ilations  showed.  This  may  be  due  to  the 
rapid  reduction  of  the  blade  load  near  the  tip  and  the  shifting  of  low 
energy  air  toward  the  trailing  face  of  the  blades  observed  for  all  test 
rvois. 


The  experimental  techniques  developed  in  the  course  Of  this  invest!* 
gat ion  can  be  used  to  obtain  consistent  and  reproducible  surveys  of  the 
internal  flow  distribution  in  rotating  cascades. 
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I.  DEDUCTION 

Considerable  progress  has  been  aehlered  in  the  past  for  the  one  and 
two^dinenslonal  methods  of  analysis  as  applied  to  radial-flow  and  axial-flow 
turbomachlnes.  However*  in  mixed-flow  machines*  the  three-dimensional  ef¬ 
fects  are  of  relatively  greater  importance*  and  further  knowledge  is  re¬ 
quired  of  the  detailed  flow  conditions*  Recent  analysei  (for  example* 
those  of  Vu  and  Staciti^*  of  inviseid*  three-dimensional  flow  Suggest 
lethods  for  the  design  of  blade  cascades.  However*  general  use  of  these 
methods  is  restricted  because  of  the  coiy>lexity  and  great  expenditoe  of 
effort  involred  in  their  application^  and  because  of  the  lack  of  experi¬ 
mental  knowledge  necessary  to  establish  the  limitations  of  the  theory. 
Iherefore*  the  objective  of  this  inrestigatien  is  to  provide  additional 
experimental  information  pertaining  to  the  flow  conditions  within  the  Im^ 
poller  and  to  evaluate  the  theory  in  terms  of  the  experimental  data.  In 
an  effort  to  siq>plenent  the  existing  experimental  techniques*  a  three- 
dimensional*  spherical  pitot-probe  and  a  pressure-transfer  device  employ¬ 
ing  mercury  seals  have  been  developed  for  measuremwt  of  the  relative  flow 
fields  within  the  rotating  impeller* 

An  experimental*  four-blade*  mixed-flow  impeller  was  constructed  in 
accordance  with  a  design  procedure  developed  from  a  bound-vortex  theory  of 
the  flow  presented  by  Spannhake^*  Axial  symmetry  was  assumed  in  the  de¬ 
sign*  and  the  Aow  was  considered  to  be  nonviscous*  incompressible*  and 
steady*  The  experimental  investigation  of  the  impeller  was  conduoted  in 


for  all  numbered  references*  see  bibliography 


MOTE:  This  publication  was  released  in  February  1956* 
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t  Qoapressor-easeade  t«8t  stmd  using  sir  as  the  fluid  medium. 

A  quaslothree-^dlmenslonal,  relaxation  solution  of  the  blade-to^blade 
▼ariatlon  in  the  flou  on  the  mean  stream  surface  of  reyolutlon  uas  cdm» 
pnted  according  to  methods  proposed  bgr  Stanlts^*^  in  order  to  investigate 
the  suitability  of  numsrlcal  methods  as  applied  to  the  more  general  curved 
blades  of  the  experimental  impeller,  the  circumferential  variation  in  the 
flow  predicted  on  the  basis  of  this  theory  is  compared  with  tbs  measured 
results* 

n.  THBDRBTICAL  AMALIS^ 

Qeperal  Theory 

The  dynamic  equation  governing  the  three-dimensional  flow  of  an  in- 
viscid  fluid  through  a  tturbonaehine,  neglecting  body  forces  > 


follows  directly  from  Newton's  second  law  of  motion*  Considering  the  den¬ 
sity  to  be  a  function  of  the  pressure  only>  the  pressure  gradient  may  be 
expressed  as  an  integral  term  by  means  of  the  relation. 


(2) 


With  reference  to  a  cylindrical  coordinate  system  fixed  to  the  impeller 


and  rotating  with  a  constant  angular  speed,  the  fluid  acceleration  is  given 


by  the  equation 


a  •  +22x7  (3) 

where  r  is  the  cylindrical  radial-coordinate,  and  Z^xw  is  the  Coriolis 
acceleration.  The  relative  accelwation  may  be  expressed  3h  tfips;  Of  the 
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relatliro  velocity  w  as 

-  dv  . 

•rel  "  ♦ 


V 


(U) 


Gonsidering  the  relative  motion  to  be  steady,  ^  -  0,  and  making  use  of 
the  vector  identity,  (W-V)  w  basic  dynamic  equation 


becoiRee 

7(1^+  ♦  2Wxw  -60*  ?  -  ir  0. 


(5) 


the  las  of  cosines  applied  to  the  triangle  forced  the  velocity  dlagrim 
ehoim  in  Pig.  1  yields  the  relation. 


w*  -  c*  ♦  (rod*-  2  r(jc^. 


(6) 


and  by  substitution  into  equation  (S),  making  use  of  the  vector  Identity 
7r*  •  2  f ,  there  is  obtained  the  equation 


It  i.  coBwnUirt  to  intpodue.  tho  «»rgT  Iwu^  H  wO  qoMttV  t  M 
cording  to  the  definitions 


C8) 


(9) 


WADC  tR  55^158 


3 


Fi|.  1  m/ocin  mmjM  vnyrnm 
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IqaatloB  (7)  lugr  be  expreiMd  as 

g7I  ♦  25X¥  -  -  0. 

the  vortlelty  of  the  absolute  aotlon  is  ralatad  to  the  ralatiire  Yortielty 
by  the  relationship 


-  X-  2uj, 

and  an  alternilte  forn  of  equation  (7a)  la 
g7I  -  wx  X  -  0. 


(10) 


(7b) 


Differentiation  of  the  quantity  I  with  respect  to  the  relative  fluid  uotlon> 
together  with  equations  (6),  (8)«  and  (9)«  leads  to  an  inportant  general 
eonsideratlont 


(gl)  • 


w  ♦  2wxu 


(11) 


For  steady  relative  uotion,  the  pressure  at  any  fixed  point  relative  to 
the  livelier  can  he  considered  constant,  and  hence,  from  equation  (11), 
the  value  of  I  for  any  fluid  particle  passing  through  the  blade  cascade 


*  A  more  general  derivation  of  948  equation  for  an  ideal  gas,  involving 
heat  conduction,  is  given  by  Wu^, 


^  -  w  X  (Vxw)  ♦  2Wxw  ■  ^  g71  ♦  77(speeiflc  enUopy), 
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will  rciwtn  constant.  In  terns  of  the  whirl  ■onentf  a  «  rc^^  the  well 
known  taler's  turbine  relation  is  obtained, 

VH  -  |  vn,  (12) 

If  the  flow  upstream  of  the  blade  cascade  has  a  uniform  distribution  of  I, 
SO  that  the  gradient  of  X  is  sero,  than,  from  equation  (Tb),  the  absolute 
vortex  lines  must  coincide  with  the  relative  streamlines,  and  the  absolute 
vorticiiy  of  the  flow  will  remain  unhanged  on  passing  throu^  the  blade 

cascade* 

B.  Infinite  Humber  of  Blades 
1.  Vortex  Theory  of.  Flow 

A  considerable  simplification  in  the  analysis  Of  the  Aow 
through  blade  cascades  is  obtained  by  considering  the  fluid  to  be  perfectly^ 
guided  by  the  blade  surfaces.  That  is  to  say,  every  relative  stream  sur* 
face  is  assumed  to  eonfon  exactly  to  the  blade  ^aqpe.  However,  this  can 
be  achieved  only  if  the  spacing  between  the  blades  is  very  small,  confining 
the  fluid  motion  essentially  to  the  curved  surface  of  the  blade.  As  ah 
Idealization,  an  infinite  number  of  blades  of  infinitesimal  thickness  will 
be  assinned#  Under  certain  conditions,  it  has  been  found  that  analyses 
based  upon  this  Idealization  will  lead  to  solutions  that  give  a  good  ap« 
proximation  of  the  mean  flow  between  blades  (Spannhake^).  ^  the  intro*^ 
duction  of  an  infinite  number  of  blades,  the  general  three-dimensional 
motion  of  the  fluid  is  reduced  to  a  two-dimensional  motion  on  the  blade 
surface  completely  described  by  two  independent  variables.  Sy  this  as- 
suaption,  the  variable  9  can  be  eliminated,  as  there  is  no  Variation  about 
the  axis  of  the  cascade,  and  the  cylindrical  coordinates  r  and  z  will 
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■ufflM  ai  tht  Indapandent  vaiiabies  for  steady  ■otlon* 

Iq  an  actual  iapeller  with  a  finite  nunber  of  blades  the  flow 
is  acted  upon  by  pressure  differences  developed  by  the  blade  surfaces. 

For  a  given  change  in  angular  mooentum  of  the  through  flow,  these  pressure 
differeneea  will  diainiah  as  the  blade  number  increases;  however,  as  the 
space  between  the  blades  gets  smaller,  the  weight  flow  acted  upon  by  the 
pressure  difference  is  reduced.  In  the  limit,  the  ratio  of  the  blade 
force  to  the  fluid  weight  is  a  finite  value,  Cofsaonly  referred  to  as  the 
"distributed  blade  force",  and  represents  the  action  of  the  blade  cascade, 
this  fictitious  force  Ic  may  be  introduced  as  an  additional  term  into  the 
basic  dynamic  equation  (1), 

I  •  -  ^  Vp  ♦  I  .  (13) 

Carrying  this  additional  force  term  through  the  steps  leading  to  equation 
(7b),  there  will  be  obtained 

gVI  ^  wxX-  ^  •  (lb) 

The  force  Ic  is  used  to  define  the  bound  vorticity  associated  with  the 
blading, 

■C  •  -  /pwx3^.  (15) 

The  vector  ^  is  normal  to  the  blade  surfeue  since  it  represents  the  pres- 
sww  force  due  to  the  blading;  and  as  w  is  tangent  to  the  blade  surface  , 
it  necessarily  follows  from  equation  (15)  that  the  bound  vortex  lines  lie 

*  This  implies,  w*lE  ■  0. 


in  the  blade  surface#  The  sub  of  the  bound  vorticity  attached  to  the 
blading:  and  f Me  vorticity  of  the  flow  is  the  total  vorticityj 

A  ■  A^  ♦  A^  • 

The  ftree  vortices  are  acted  upon  only  by  the  gravitational  fields  they 
follow  the  Helmholta  laws,  aid  have  no  relationship  with  the  force 
Equation  (lU)  nay  then  be  expressed  as 

gVI  -  wxA^  ^  wxAj  “  ^ 

ft'OB  the  definition  of  the  bound  Vortices,  equation  (15) »  it  follows  that 
g7I  -  wxA^  -  0  .  (18) 


For  the  sane  reasons  as  discussed  ijnnediately  follcwing  equation  (11)  *  the 
fluid  particles  that  have  a  uniform  distribution  of  I  in  a  region  upstreaiii 
of  the  blade  cascade  will  move  along  relative  streamlines  which  satisfy 
the  relation 


Equation  (19)  implies  that  either  the  free  vorticity  is  zero,  or  the  free 
vortex  lines  coincide  with  the  relative  streamlines. 

2,  Special  Case  of  Zero  Ring  Vertex  (  A^  p  0) 

The  impeller  used  in  the  experimental  investigation  was  de«^ 
signed  in  accordance  with  the  preceding  theory  for  an  inviscid  fluid  pass^ 
Ing  through  a  cascade  consisting  of  an  infinite  number  of  blades.  A  fur^ 
ther  siigjlification  in  the  4esign  procedure  is  introdneed  by  tnposijig  the 
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condition  that  the  tangential  component  of  the  vorticlty  be  zero.  In  terme 
of  the  cylindrleal  coordinate  system,  the  tangential  con^onent  of  the 


vorticity  is 


dc  be 
_ r  _ z 

sr  "sr 


0  . 


(20) 


On  the  basis  of  the  continuity  equation,  it  is  necessary  that  the  radial 
and  axial  components  of  both  the  absolute  and  relative  fluid  vslocities  at 
any  point  be  identic^,  that  is. 


c  «  V.  , 
r  r  * 


and  hence. 


bw, 

T 

os  “  dr 


«  0 


(21) 


The  equation  of  the  blade  surface  may  be  represented  by  the  equation 

F(r,z)  ♦  0  -  0  .  (22) 

The  direct’iou  numbers  of  a  normal  n  to  the  blade  surface  are  (^,  #)* 

The  free  vortex  lines  must  lie  in  the  blade  surface  because,  they  coincide 
with  the  relative  streamlines.  Further,  since  the  bound  vortex  lines  lie 
in  the  blade  surface,  it  follows  that  the  total  absolute  vortex  lines  must 
also  lie  in  the  blade  surface  as  they  are  deterslned  by  the  vmetor  sum  of 
the  fbee  and  bound  vortices*  This  implies  that 


and 


n  ^  X  ^ 
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the  infinite  blade  restriction  Implies  axial  eyinmetry  so  that  differentia¬ 
tion  with  respect  to  0  is  zero,  introducing  the  whirl  npmenty  ■  ■  re^* 
and  for  ■  0,  equation  (23)  sliqilifies  to 


6F  dn  dF  6bi 

^  ”  Sr 


0  , 


m 


The  left  side  of  this  equation  is  the  Jacobian  determinant  of  the  functions 
F(r>z)  and  m(r>s),  and  the  ranishing  of  this  determinant  iinplies  the  func¬ 
tional  dependence  of  F  and  m.  This  means  that 


n  •  m(F) 

^  »(9),  (25) 

andy  consequentljTy  the  Intersection  of  the  blade  surface  with  radial  planes 
(9  m  constant)  are  lines  of  constant  whirl-moment. 

3.  jBlade  Design  for  Zero  Ring-yortex 

a*  Petermlnation  of  the  Meridional  Streamline  Pattern 

As  a  preliminary  consideration^  suppose  there  are  no  blades 
present  between  the  hub  and  shroud  surfaces*  Irrotatlonal  flow  (i*e*  flow 
haring  aero  Torticity)  passing  through  this  annular  space  will  satisfy  the 
condition 
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I 


(26) 


sr 


'5r 


This  diff«r«ntlal  equation  may  be  used  as  the  basia  for  the  meridional 
ttreamllne  pattern  of  the  free  flow  through  the  annular  apaee  without 
blading.  If  bladesi  or  the  equlTalent  bound-vortex  surfaces j  are  Intro¬ 
duced  such  that  the  tangential  component  of  vortlclty  Is  iero»  then  equa¬ 
tion  (21)  will  be  valid.  This  means  that  equation  (26),  which  applies  to 
the  free  flow  without  blading,  will  remain  unchanged  by  the  presence  of 
blades  deiigned  so  that  Is  identically  zero. 

The  first  step  In  the  blade  design  Is  the  determination  of 
the  meridional  flow  pattern  for  free  Irrotatlonal  flow  between  specified 
contours  of  the  hub  and  shroud  surfaces.  The  natural  coordinate  system 
(s,n,0)  is  convenient  to  use  for  this  determination*  The  flow  Is  irroti- 


tiooal. 


i>c_  be 


5r 


(27) 


and  the  convective  component  of  the  acceleration  in  the  direction  normal 


to  the  streamline  Is 


c 


m  OS 


(28) 


and,  hence,  the  Irrotatlonal  oondition  may  be  ejqpreseed  by 


(29) 


This  differential  equ*^loh  may  be  integrated  iLong  an  orthogMial  (llg*  2), 


11 
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In 


(30) 


The  velocity  distribution  along  an  orthogonal,  given  by  equation  (30)* 
satisfies  the  Izrotatlonal  flow  condition.  If  the  meridional  streamline 
pattern  is  constructed  to  satisfy  this  condition,  together  with  a  contin¬ 
uity  requirement,  the  streamline  pattern  will  be  uniquely  determined*  The 


continuity  condition  requires  that  the  same  flow  rate  cross  any  two  sec¬ 


tions  between  adjacent  stream  surfaces.  The  volume  rate  of  flow  passing 
between  the  hub  surface  (n  ■  0),  aid  any  other  stream  surface  is 


and,  in  dimensionless  form. 


(31) 


Equations  C30)  and  (31)  can  be  satisfied  by  a  graphical 
method  of  conqiutation.  For  the  specified  hub  and  shroud  surface  contours 
of  Fig.  2,  a  number  of  estimated  orthogonals  are  sketched  into  the  mer¬ 
idional  plane*  For  each  of  the  orthogonals,  the  velocity  distribution 
given  by  equation  (30)  is  solved  ty  graphical  integration.  This  is  done 
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most  easily  by  estimating  the  distribution  of  Lti/t  along  the  orthogonal, 
and  plotting  this  distribution  against  n/An.  The  area  under  this  curre  up 
to  any  value  of  n/^n,  represents  the  dimensionless  velocity  at  this 

point;  and  from  the  area  Values,  the  distribution  shown  In  Fig*  3  is  ob* 
tained*  The  dimensionless  discharge,  also  shown  in  Fig*  3*  is  computed 
from  equation  (31)  by  a  graphical  integration.  From  this  discharge  plot, 
the  n-coordinates  of  the  streamlines  which  equally  partition  the  flow 
along  each  orthogonal  are  determined.  Using  the  appropriate  radius  and 
center  of  curvature  at  each  of  these  n-coordinates,  the  arcs  of  the  mer¬ 
idional  streamlines  azs  drawn  in  with  a  compass*  If  the  arcs  blend  into 
a  smooth  curvei  the  correct  meridional  streamlines  have  been  obtained^ 
satisfying  both  the  continuity  and  irrOtational  flow  conditions*  If  the 
arcs  do  not  blend  together,  the  original  estimates  of  the  orthogonals  must 
be  modified  and  the  conputational  procedure  repeated*  The  dimensionless 
leyout  of  the  meridional  streamline  pattern  shown  in  Fig.  2  was  obtained 
after  the  second  confutation*  The  pertinent  numerical  velues  for  this 

streamline  pattern  are  presented  in  Table  X*  The  distribution  of  the  di- 

®m 

mensionless  velocity  -  -»«  along  each  orthogonal,  and  along  each  meridional 
streamline^  is  shown  In  Figs.  U  and  5  respectively* 
b*  Petermination  of  the  Blade  Shape 

The  design  of  the  blade  surface,  considered  as  a  bound 
vortex-sheet,  is  such  that  at  the  design  operating  condition  the  ring- 
component  of  the  bound  vortlcity  is  zero*  The  design  operating  condition 
is  dependent  upon  the  desired  discharge  for  a  given  speed  and  impeller 
size,  and  is  fixed  by  a  suitable  choice  of  the  design  parameter,  Q^iD^* 
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Zt  n*  76  To  1.0 

Orthogonal  Coordinato  ^ 


Pig.  3  DISCHARGE  AND  MERIDIONAL  VELOCITY  DISTRIBUTION 
ALONG  ORTHDGONAL. 


WifflC  TR  55-158 


15 


Dlinensloniess  Discharge  > 
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The  blade  design  le  based  \xpon  a  typical  Talus  for  this  type  of  mixed-flow 
Impeller  of  ^0^  ■  0.060.  The  blade  cascade  satisfying  these  conditions 


will  not  disturb  the  meridional  streamline  pattern  already  determined  for 


the  free  flow. 

Equation  (25)  expresses  the  fact  that  the  Intersection  line 
of  the  blade  surface  with  any  radial  plane  Is  a  line  of  Constant  whlrlnnoment. 
The  whirl-moment  Is  a  function  of  0,  and  can  be  arbitrarily  Chosen;  how- 
everj  seTsral  trials  may  be  necessary  to  produce  a  desirable  blade  shape. 

For  example*  the  total  length  of  the  blade  may  turn  out  to  be  Tery  different 
along  the  hub  and  the  shroud  surfaces.  How  to  change  the  whirl-moment  for 
a  desired  modification  of  the  blade  shape  is  easily  understood*  once  one  has 
gone  through  the  design  procedure.  The  first  choice  of  whirl-moment  for  the 
blade  design  of  the  test  iBg>eller  led  to  a  double- currature  of  the  blade 
along  the  shroud  surface.  The  whirl-moment  was  adjusted  so  that  this 
double-curvature  was  eliminated*  and  the  final  choice  of  whirl-moment  dis¬ 


tribution  is  shown  In  Fig.  6. 

The  design  of  the  acperimental  blade  cascade  was  acconqpllshed 
by  a  combined  numerical  and  graphical  procedure.  The  differential  equation 
of  the  relative  streamlines  on  a  conical  development  is 


4  -  r  H  -  cot  p , 
n 


(32) 


and  may  be  approximated  by  a  finite-difference  equation  expressed  In  di¬ 
mensionless  form* 


60  «  ^  Qotpas. 


(33) 
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XroB  this  equation,  angular  Increments  were  computed  for  constant  Incre¬ 
mental  arc  lengths,  Ss  ■  0.028,  starting  from  the  blade  exit  edge  chosen 
idong  orthogonal  V  as  shown  In  Fig.  2.  The  corresponding  values  of  R  were 
measured  directly  from  this  figure,  and  the  values  el  eotp  were  determined 
from  the  velocity  diagram  at  each  point.  It  Is  apparent  from  the  velocity 
diagram  shown  In  Fig.  1,  that  the  tangential  and  meridional  relative  vel¬ 
ocity  components  are  given  by 

^  -  °t  » 

•m-^'m  • 

Introducing  the  whlrl-momentf  m  >  rc^,  the  value  of  cot^  IS  determined 
from  the  equation 

u-f 

cotp  •  '  '  , 

m 

and,  In  dimensionless  form. 


cotp 


Q/D 


(3U) 


The  dimensionless  meridional  velocity  conq>onent  Is  taken  from  Fig.  and 
the  whirl-moment  from  Fig.  6.  The  details  of  the  design  procedure  are 
progressively  shown  In  Fig.  7,  and  the  numerical  values  for  the  computa¬ 
tion  are  given  In  Table  II.  The  constant  Incremental  arc  lengths  Ss 
along  each  meridional  streamline  are  shown  In  Fig.  7a.  The  corresponding 
co^juted  values  of  are  used  to  determine  the  streamlines,  shown  In  the 
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(«)  Cons  true  tlon  of  feints  IMsplacsO  from  Blsdo  Outlot 

Kdgs  A  Constant  Arc-Distance  Alone  Meridional  Stroanllnas 

rt(.  7  GRAPHICAl  BUDI  DCSICN  niOCBOURS 


WAPe  1!R 


Angls  betuten  leiding  and 
trailing  blade  surfaces 


(b)  Plan  View  of  Hsridlonsl  Strsanainss  and  Radial  Planes 
Pig.  7  (Continued)  QRAPKICAL  BLADE  DESIGN  PROCEDURl 
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tterldloml  ■trotmllnel 


(c)  Radially  ?ut  Blade  Intersection  lines 
ri*.  7  (Concluded)  ORAPHICAL  BLAilS  tBSICN  BROCBWRE 


plan  view  of  ?lg.  7b.  Starting  from  the  blade  outlet  edge,  radial  planes 
«re  drawn  In  the  plan  view  at  ^  degree  Intervals  up  to  the  inlet  blade 
edge,  spanning  a  total  cngular  displacement  of  $0  degreea.  The  intersect 
tion  points  of  the  meridional  streamlines  with  these  radial  planes  are 


plotted  in  Fig.  7c,  and  these  points  determine  the  Intersection  lines  of 
the  blade  surface  with  the  radial  planes.  As  previously  described,  these 


intersection  lines  are  lines  of  constant  whlrl>mo!(ient,  and,  consequently. 
Indicate  lines  of  constant  energy  input  to  the  fluid  from  the  blade  Cas-^ 
cade,  the  eyllndrlcal  coordinates  of  the  blade  surface  are  measiured  di¬ 
rectly  f**om  these  intersection  lines  and  are  presented  in  Table  III. 
Tesplates  designed  according  to  the  data  of  Table  III  were  used  to  guide 
the  construction  of  the  pattern  for  the  impeller  easting. 

C.  Finite  Number  of  Blades 
!•  General  Considerations 

The  design  of  the  e}q)erlmental  impeller  is  based  upon  flow 
conditions  through  a  cascade  with  an  infinite  number  of  blades.  Ruden^ 
has  computed  the  potential  through-flow  for  a  cascade  having  a  finite 
number  of  blades  by  considering  the  velocity  at  every  point  to  be  the  sum 
of  the  average  value  with  respect  to  0  plus  a  perturbation  velocity.  He 


showed  that  if  perturbations  of  an  order  higher  than  the  first  are 


neglected,  the  average  meridional  values  with  respect  to  0  «re  exactly 
equal  to  the  values  computed  on  the  basis  of  an  Infinite  number  of  blades. 
In  general  practice,  this  condition  is  considered  to  be  satisfied  If  the 
blades  are  not  spaced  too  widely  apart,  the  fluid  is  of  low  viscosity, 
and  the  operating  conditions  are  near  the  design  point.  One  of  the  pb- 
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objeetlves  of  the  experimental  phase  of  this  program  is  to  determine 
the  esctent  to  vhleh  the  flow  deviatee  from  the  theoretically  predicted 
flow  for  an  infinite  number  of  blades*  Assuming  a  uniform  flo-r 
distribution  from  hub  to  shroud,  Stanitz^*^  has  computed  the  blade-to*^blade 
variation  in  the  flow  through  a  conical,  mixed-floif  compressor  by  the  re¬ 
laxation  me^od.  He  has  also  conQ)Uted  the  flow  on  a  general  mean  stream 
surface  of  revolution  for  a  compressor  having  strai^t  radial  blades. 
Following  the  methods  of  Stanitz,  the  theoretical  variation  of  the  flow 
between  the  blades  of  the  experimental  impeller  has  been  computed  at  the 
design  operating  condition  Q/(aD^  •  0.060,  for  comparison  with  experiments 
and  evaluation  of  the  underlying  potential  flow  assumptions. 

Because  of  the  small  blade  number  (four)  of  the  expeFimental 
lii9)eller,  the  maximum  deviation  from  the  mean  flew  conditions  are  expected 
to  occur  in  the  tangential  direction,  the  spacing  between  the  hub  and 
shroud  surfaces  is  relatively  small  compared  to  the  spacing  between  the 
blades;  thus,  it  is  assumed  that  the  deviation  of  the  meridional  streamline 
pattern  from  the  axlal-synaetry  solution  is  small*  Further,  at  the  design 
condition,  the  tangential  component  of  the  bound  vortex-sheets  representing 
the  blade  surface  has  been  assumed  to  be  zero,  so  that  from  the  infinite- 
blade  considerations  the  meridional  streamline  pattern  should  remain  un- 
distvebed.  For  these  reasons,  the  flow  is  assumed  to  proceed  along  the 
surfaces  of  revolution  generated  by  the  meridional  streamlines*  This  as- 
svpption  does  not  allow  a  true,  three-dimensional  motion,  since  the  fluid 
particles  are  thus  hypothetically  confined  to  predetermined  surfaces  of 
revolution*  The  quasi-three-dimenslonal  solution  so  obtained,  however, 
should  give  a  more  realistic  picture  of  the  flow  conditions  within  the 
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lipelltr  than  would  thi  «xi*l-8ynnotr7  solution. 

2*  Egustlon  of  Fluid  Motion  on  the  Moan  Msfldleftsl  Strwsm  Surfses 
a.  Al?soluta  Irrotatlonal  Motion 

Tha  ganaral  differential  equation  of  the  fluid  motion  Is 
derived  under  the  assvuqjtlon  that  the  flow  is  confined  within  a  ehannel 
formed  by  two  stream  surfaces  of  revolution  enclosing  the  stream  surface 
generated  by  the  mean  meridional  streamline  as  shewn  in  Fig.  6.  the  flow 
Is  assumed  to  be  Irrotational  upstream  of  the  blade  cascade «  and  according 
to  tha  discussion  immediately  following  equation  (12),  the  absolute  velocity 


of  the  fluid  will  remain  Irrotatlonal  upon  passing  through  the  Iflqjeller. 

In  particular,  the  vortlcity  component  normal  to  the  stream  surface  Is  aero. 


Implying  the  circulation  on  any  closed  path  on  the  stream  surface  between 


the  blades  is  tero  (Stokes'  theorem).  Considering  the  path  shown  In  Fig. 
the  following  equation  is  obtained, 

<«>«]  • 

Since  sinS  «  dR/dS,  this  equation  may  be  simplified, 

bw 

[rC2R  ♦  w^)]  -1^3^  j 
and,  smpandlng  this  equation  gives 

^  dW.  dW. 

liR  sinS  ♦  sinS  ♦  R  slno*^  ^  3^  •  0  .  (36) 

b.  Continuity  Egustlon 

The  application  of  continuity  oondltlona  to  tha  alamantal 


WADC  TR  55-158 


27 


Circulation  about  path 

dT-  iU(c,ii(ie)dH-  ;^e_ds):<ie 


Fig.  9  GIRCmjlTION  ANB  vm/JCITT  DlAOiAH  ON  MEAN  STREAM  SR^CE 


ngioQ  thowa  In  Pig.  10  ylnlds  th«  following  oquatlon. 

0  «  ^  (V^ER  dO)dR  ^  (V^B  dS}dO  ; 

«ndf  itnco  dR/dS  ■  ainS  ,  thoro  is  obtalnad  the  aquation 


at  <«.»>  ♦  ^  <*'t 


■  0  . 


A  straan  function  <//  .  satisfying  equation  (37) >  may  be  deflnad  such  that 


B 

sind  * 


(38a) 


4jf  ^  •  (38b) 

Coordlnatas  (^.iji).  introduced  bgr  Stanlti^.  on  tha  maan  stream  surface  are 
deflnad  such  that 


1 

sz  * 


mhm  • 


(39a) 

(39b) 


Since  the  Cauchy-^Rlemaim  equations  are  satisfied,  the  transformation  of  co^ 
ordinates  (R.  0)  fron  the  mean  stream  avrfaea  to  the  (C  fV)  plane  la  eonr 
formal}  that  la. 


d(  idV 
5S"55^  ' 

1^  ^ 

In  the  transformed  plane,  equation  (36)  becomes 


(UOa) 

(UOb) 
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1  dW.  .  6W^ 

l4H-*in6  ♦  8in8+ ■  0, 

•Dd  tht  •quations  for  th«  stream  funetlofi  become* 


(Ui) 


(Use) 


^  BRAO, 


e.  Qeneral  Eqtaatlon 

Combining  the  etream-f unction  equations  (k2)>  and  equation 
(UI),  the  equation  given  by  Stanlts  Is  obtalnedf 


4^  *  4^  ■  Tf  ■  BE*(68)*  slnS  .  0  . 

Using  the  notations  shown  In  ^g.  11»  equation  (U3)  nay  be 
flnlteHllfferenee  form, 

^2*  ^3*  ^  ^0^  ^  ^^3^  ■  **^*^^0^  “ 


where 


g(^)  -  u[bR^(AO)^  slnS]  . 


m 

in 

(UU) 

(U5) 

(U6) 


The  values  of  ^  as  a  function  of  R  are  determined  from  a  numerical  evalua* 
tlon  of  the  Integral 


XV 

1  f  1 

"A^y  Fsi 


^lade  exit 


a?) 


The  plot  of  the  mean  meridional  streamline  and  the  channel  height  6  are 
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niiminiiwiini 


shoim  In  Fig.  12.  The  eohe  angle  B  ie  determined  from  the  slope  of  the 
mesa  merldionel  etreamline  with  respect  to  the  ucis  of  rotstlonf 

6  -  arc  ten  ^  ,  (1*8) 

and  the  values  of  sin 8  are  plotted  in  Fig*  13*  Ths  values  of  C  obtained 
by  the  numerical  integration  of  equation  (1*7)  are  shown  in  Fig*  ly,  and 
the  function  InB,  together  with  its  derivative,  is  shown  in  Fig.  11* •  the 
0-eoordinates  of  the  blade  profile  on  the  mean  Stream  surface  are  shoim  in 
Fig*  12,  and  the  coxTesponding  ^-coordinates  are  plotted  in  Fig*  13*  the 
numerical  values  are  presented  in  Table  TV* 
d*  flpundajy  Conditions 

The  assumed  boundary  conditions  that  make  the  problme 
completely  determinate  cm  be  summarised  by  the  following  four  conditions* 
1*  The  specified  inlet  flow  is  uniform  far  igMtream* 

2*  The  Kutta-Joukowsky  condition  ippliss  at  the  blade 
exit  tip* 

3*  Periodic  flow  conditions  prevail  about  the  impeller 
axis,  with  a  period  equal  to  the  blade  number* 

1**  The  flow  is  uniform  far  dounstream  of  the  blade 
cascade. 

The  flow  conditions  in  the  draft-tube  upetream  of  the 
blade  cascade  are  considered  uniform  at  C  «  1*10*  Velocity  traverses 

made  during  the  test  runs  shoved  this  assumption  to  be  essentially  correct* 
For  this  value  of  the  radius  R  is  0*2^0,  and  the  channel  height  B  is 
2*00*  The  meridional  component  of  the  upstream  velocity  is  given  by  the 


WADC  IR  55-158  32 


Fig,  12  MEAN  MERIDION/IL  StRE^OlMl  CgOM31NATES  AND  CHANNEL  HEIGHT 
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1 

R 

B 

-  1.10 

0.250  ' 

2.00 

1.00 

0.253 

1.95 

-  0,90 

0.255 

1.89 

-  0,80 

0.258 

1.85 

^  0.70 

0.265 

1.78 

-  0,60 

0.27li 

1.71 

-  0-.50 

0.287 

1.62 

-  o.Uo 

0.306 

1.52 

-  0,30 

0.331 

1.39 

<-  0,20 

0,367 

1.26 

0.10 

O.I4I6 

1.12 

*  0*05 

0,i4lt6 

1.06 

0 

O.I477 

1.00 

0.05 

0.511» 

O.9U 

0.10 

0.551 

0.90 

0.15 

0.596 

0.87 

0.20 

0.6141 

0.85 

0.25 

0.691 

O.8I4 

0.30 

0,7l45 

0.814 

0 

0.065 

0.093 

0.125 

0,187 

0.265 

0.365 

0.1490 

0.620 

0.7li8 

0.860 

0.907 

0,9142 

0.968 

0.982 

0.995 

1.000 

1.000 

1.000 
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The  span  between  blades  is  AT}*  and  consequently,  0.2276*  For 
convenience,  a  stream-function  may  be  defined  such  that 

>1/  .  100  .  (SO) 

A\l/ 

The  value  of  ^  on  the  leading  face  of  the  blade  profile  is  made  equal  to 
aero,  and  on  the  trailing  face  a  value  of  100  is  chosen.  The  specified 
values  of  the  derivatives  of  the  stream-function  on  the  v^stream  boundary 
are  shown  in  Fig.  15« 

The  Kutta-Joukowsky  condition  requires  that  the  relative 
flow  leaves  the  blade  tip  tangent  to  the  blade  surface.  The  best  numerical 
procedure  to  follow  is  to  satisfy  all  boundary  conditions  with  the  excep¬ 
tion  of  the  Kutta  condition  in  the  first  relaxation  solution,  and  then  ad¬ 
just  the  downstream  whirl-moment  value  until  the  tangency  condition  is  sat¬ 
isfied.  The  adjustment  of  the  downstream  whirl-moment  value  is  determined 
by  a  separate  relaxation  solution  involving  fixed  zero-values  of  the  stream- 
function  on  the  blades  and  a  unit  change  in  the  upstream  whirl-moment  value. 
The  corresponding  changes  in  the  stream-function  are  determined  to  a  first 
s^proximatlon  by  the  relaxation  solution  of  the  Laplace  equation, 

0,  (51) 

and,  in  accordance  with  these  values,  the  stream-function  is  adjusted  to 
satisfy  the  Kutta  condition.  A  final  coiqputation  must  then  be  made  to 
satisfy  the  more  accurate  differential  equation  (U3),  The  details  of 
this  computation  are  presented  in  Appendix  C. 

Prom  symmetry  considerations,  the  flow  is  cyclic  in  the 
-direction  with  a  period  equal  to  the  blade  spacing,  A‘7»  1.  Periodic 

3rr 
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are  constant  on  downstream  boundary 


0  (leading  100  (trailing 

blade  face)  blade  face) 


L-rj^l.O 


^  Periodic 
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^  Upstream  boundary  ^.t. 


Pig.  15  BOIIKDAPJ  CONDITIONS  FOR  THE  FLOW  FIELD  ON  THE  MEAN  STREAM 
SURFACE. 
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boundary  lines  are  eketched  in  as  shown  In  Fig.  1$,  and  for  any  given 
value  of  the  value  of  the  stream-function  on  the  right-hand  boundary 
line  must  equal  tiie  value  on  the  left-hand  boundary  line  plus  100. 

The  flow  conditions  are  Considered  to  become  uniform 
downstream  in  the  vaheless  diffuser^  so  that  the  derivatives  of  the  stream^ 
function,  shown  in  Fig.  15,  are  constant  on  the  downstream  boundary. 

The  flow  conditions  around  the  boundaries  enclosing  the 
flow  region  have  now  been  completely  described,  and  the  problem  can  be 
solved  by  relaxation  methods.  A  discussion  of  the  results  of  this  ooi^u- 
tatlon  at  the  design  condition  follows. 

3.  Results  of  the  Jbelaxatipn  Solution 

After  the  residuals  oCcuring  in  the  relaxation  solution 
have  been  eliminated,  the  values  of  the  stream  function  at  the  nodes  of 
the  relaxation  grid  are  used  to  ctetermlne  the  relative  streamlines.  A 
streamline  is  a  line  drawn  smoothly  through  coanon  values  of  the  stream 
function  as  is  shown  in  Fig.  16.  For  any  streamline,  the  numerical  value 
of  ^  indicates  the  quantity  of  flow,  passing  between  the  contour  line  of 
the  leading  blade  face  and  the  specified  streamline,  as  a  percentage  of 
the  total  flow  passing  between  the  two  blades.  From  the  distribution  of 
the  stream-function  along  the  periphery  at  the  blade  exit-tip,  the  mer¬ 
idional  velocity  distribution  is  coa^uted  by  means  of  equation  (U2b). 

The  velocity  distribution  con5)uted  on  this  basis  is  coBqpared,in  Chapter 
IV,  with  the  distribution  measured  at  various  operating  conditions  of  the 
test  li^Mller. 
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m.EIPERIMBWrAL  INVESTIQATIOW 

A*  Description  of  the  Test  Stand 

A  schematic  diagram  of  the  test  stand  used  in  this  investigation 
is  shown  in  Fig.  I?,  and  the  conplete  test  stand  with  ly^e  dynamometer 
drive  is  shown  in  the  photograph  of  Fig.  18.  The  dynamometer  unit  has  a 
20  horsepower  rating  and  a  speed  range  from  1^00  to  6000  rpm.  A  large^ 
box-type  filter  is  provided  upstream  of  the  inlet  nozzle  to  prevent  damage 
and  clogging  of  the  instruments.  The  filter  consists  of  flannel  cloth  sup¬ 
ported  by  wire  mesh  screening^  and  has  static  pressure  taps  at  the  walls 
of  the  box  to  measure  the  pressure  drop  across  the  filter.  A  standard 
A.S.K.E.  long-radius  nozzle  is  installed  at  the  inlet  of  the  draft  tube  to 
measure  the  rate  of  flow,  and  to  provide  smooth  entry  of  the  air.  Six 
straightening  vanes  are  mounted  in  the  nozzle  (Fig.  19),  and  extend  from 
the  inlet  face  to  the  end  of  the  long  radius » 

The  impeller  is  cast  of  an  aluminum  alloy  and  has  a  maximum  blade 
diameter  of  1U.Q6  inch  (see  Figs.  20  and  21).  The  details  Of  the  Impeller 
design  have  been  described  in  Section  II.  An  aluminum  disc  is  mounted  on 
the  bach  of  the  impeller  to  s\;pport  the  pitot  probe,  which  can  be  located 
in  six  different  peripheral  positions  relative  to  the  impeller.  An  ad¬ 
justable  dimny  pitot-probe  is  provided  diametrically  opposite  the  actual 
probe  to  minimize  the  dynamic  unbalance  of  the  unit.  The  impeller  and 
the  positioning  disc  were  each  dynamically  balanced  so  that  the  unit  remains 
in  dynamic  balance  even  though  the  location  of  the  positioning  disc  is 
changed  relative  to  ths  Impeller.  The  ring  arringement  shown  on  the  back 
of  the  positioning  disc  is  to  support  the  stainless*steel  tubing  which 


kl 


Fig,  17  SCHEMRTIC  DULGEAM  OF  THE  TEST  STAND 
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Pig.  18  TEST  STAN©  AND  DYNAMOMETER  UNIT 


Jig.  19  INLET  NOZZLE  AND  STRAIGOTENING  VANES 


“t 


leads  from  the  pitot-probe  to  the  shaft,  against  centrifugal  force  effects. 
The  tubing  is  connected  to  longltudlr^l  holes  which  are  drilled  In  the 
Shaft  and  lead  to  the  rotary-aeaL  pressure  take-offs. 

At  the  exit  of  the  Impeller,  a  vaneless  diffuser  Is  provided  whlGh 
consists  of  two  parallel  plates  of  large  diameter.  Qulde  vanes  are  omitted 
so  that  the  relative  flow  in  the  Impeller  will  be  undisturbed  and  steady. 
The  flow  leaving  the  diffuser,  passes  Into  a  discharge  drum  which  is 
equipped  with  adjustable  louvers  to  control  the  Impeller  discharge  pressuire 
and  the  total  rate  of  flow  passing  through  the  unit. 

6.  ins  truae  n tat ion 

1.  Inlet  Gondltlons 

The  total  volume  rate  of  flow  entering  the  draft  tube  Is 
measvired  with  a  standard  A.S.M.E.  long-radius  nozzle.  The  static  pressure 
taps  are  in  the  parallel  throat  portion  at  the  recommended  distance  from 
the  inlet  face.  Values  of  the  nozzle  coefficient  for  a  range  of  throat 
Reynolds  nunfcers  are  provided  in  the  A.S.M.E.  Fluid  Meters  Report. 

Conditions  downstream  of  the  straightening  vanes  in  the  nozzle 
are  measured  by  means  of  a  pitot-traverse  taken  with  a  spherical  probe. 
Pre-rotation  of  the  flow  could  not  be  detected,  either  with  the  pitot- 
traverse,  or  by  visual  observation  using  small  flags  inserted  into  the 
flow  stream. 

2.  Outlet  Conditions 

a.  Preliminary  ConaideratjonB 

The  absolute  velocity  of  the  fluid  leaving  the  impeller  is 
unsteady,  and  it  varies  periodically  with  a  frequency  given  by  the  product 
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of  the  blade  number  and  the  impeller  rotative  speed.  Because  of  the 
errors  produced  by  pressure  waves  and  high  fluid  resistance  in  the  tubing 
connecting  the  probe  and  the  pressure  pickup,  accurate  measurement  of  the 
absolute  velocity  (and  pressure)  distribution  by  means  of  stationary  total 
pressure  probes  is  considered  impractical.  There  are  pressure*^sensitive 
electronic  instruments  commercially  available  which  are  highly  responsive, 
but  an  accurate  measurement  of  the  direction  of  the  unsteady  absolute  flow 
is  considered  unfeasible  using  any  type  of  total  pressure  probe  presently 
available j  extreme  sensitivity  would  be  required  of  the  instrument,  par's- 
ticularly  at  low  flows,  and  very  small  size  would  be  necessary  to  eliminate 
disturbances.  Studies  of  incompressible  liquid  flow  In  pumps  have  been 
successfully  made  using  an  intercepting  valve,  rather  than  an  instrument 
recording  a  continuous  pressure  measiarement  with  time.  In  this  method, 
pressure  impulses  of  very  short  duration  are  intercepted  and  measured  at 
intervals  synchronized  with  the  impeller  rotation;  and,  by  varying  the 
phase  of  interception,  a  series  of  pressure  impulses  can  be  plotted  as  a 
continuous  curve.  At  the  greater  impeller  speeds,  serious  objections  to 
this  method  are  raised  because  there  is  hardly  any  means  of  actuating 
the  intercepting  valve  fast  enough  to  transmit  pressure  impulses  within 
small  phase  angles. 

A  possible  alternate  method  of  measuring  the  unsteac^ 
absolute  velocity  is  by  means  of  the  hotfwire  anemometer.  As  a  prelimin¬ 
ary  study,  small  hot-wire  probes  were  made  with  0,0U5  inch  lengths  of 
0,0003  inch  diameter  tungsten  wire  and  were  calibrated  in  a  steady  flow 
stream.  Difficulty  was  encountered  because  of  the  rather  low  sensitivity 
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of  the  hot-wlrej  when  used  as  a  constant-current  instrument,  at  air  vel¬ 
ocities  over  IhO  fps.  Greater  sensitivity  could  be  obtained  by  increasing 
the  current,  but  at  the  risk  of  burning  out  the  wire;  a  100  mill lamp  cur¬ 
rent  was  considered  to  be  the  maximum  safe  operating  limit.  It  was  decided 
that  this  instrument  was  not  practical  to  use  for  accurate  quantitative 
measurements  in  the  range  of  air  velocities  encountered,  A  further  objec¬ 
tion  to  the  hot-wire  probe,  was  the  significant  change  in  calibration 
caused  by  deposits  on  the  wire  from  small  traces  of  oil  vapor  in  the  air. 
the  hot-wire  probe  and  electronic  circuits  employed  were  of  the  type  used 
hy  Pearson^^*^^. 

Because  of  the  difficulties  involved  in  measurement  of  the 
unsteady  absolute  flow,  it  was  decided  to  measure  the  steady  relative  flow 
by  means  Of  a  pitot- type  probe  mounted  on  the  impeller,  the  two  chief  dif¬ 
ficulties  involved  with  this  type  of  probe  were:  (l)  adjustment  of  the 
probe  while  in  motion,  and  (2)  transference  of  the  pressure  from  the  ro¬ 
tating  probe  to  stationary  manometers,  A  discussion  of  these  difficulties 
and  the  manner  in  which  they  were  overcome  is  Included  in  the  following 
paragraphs. 

b.  the  three-Bimensional  Spherical  PltotrProbe 

Among  the  several  types  of  pitot-probes  commonly  used  for 
measurement  of  flow  direction,  as  well  as  flow  velocity  and  static  pres¬ 
sure,  are  the  claw- type  probe  and  the  spherical  probe,  the  most  direct 
method  of  obtaining  flow  direction  with  the  spherical  probe  is  to  adjust 
the  probe  along  two  perpendicular  planes  until  its  axis  points  in  the 
direction  of  flow,  this  direction  is  then  measured  by  protractor  scales 
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on  the  two  planes*  Coincidence  of  the  probe  axle  with  the  flow  direction 
is  obtained  by  balancing  the  pressure  readings  on  synmetrlcally  positioned 
orifice  openings.  An  alternative  method  of  measuring  flow  direction  oon- 
elsts  of  yawing  the  probe  until  the  resultant  velocity  vector  Is  In  one 
known  plane,  and  then  using  two-dimensional  calibration  curves  for  the  de¬ 
termination  of  the  flow  direction  In  the  other  plans.  Both  of  these  me¬ 
thods  require  adjustment  of  the  probe  during  operation,  which  Involves 
almost  Insurmountable  mechanical  difficulties,  If  accurate  positioning  Is 
to  be  obtained.  To  avoid  this  difficulty,  a  three-dimensional  calibration 
scheme  was  devised  so  that  any  velocity  direction  within  an  incident  cone- 
angle  of  25  degrees  can  be  determined  by  measuring  the  pressure  readings 
at  five  individual  holes  on  the  sphere  surface*  The  details  of  this  cal¬ 
ibration  are  Included  In  Appendix  B. 

The  probe  consists  essentially  of  a  tapered  shaft  with  a 
0*190-lnch-diameter  spherical  senslng-head.  The  shaft  Is  welded  at  a 
right-angle  to  a  crosspiece  mounted  between  the  hub  and  shroud  surfaces 
of  the  blade  cascade.  Five  pressure  holes  are  on  the  spherical  surface} 
the  center  hole  Is  at  the  end  of  the  shaft  axis,  and  each  of  the  other 
four  holes  are  symmetrically  spaced  on  a  sphere  radius  forming  an  angle 
of  UO  degrees  with  the  probe  axis*  Detailed  drawings  of  the  probe  are 
shown  In  Appendix  B.  Two  probes  were  made,  one  having  orifice  holes  of 
0.011  inch  diameter,  and  the  other  having  holes  of  0*016  Inch  diameter* 

The  Individual  calibration  curves  for  the  two  different  probes  are  com¬ 
pared  In  Appendix  B.  The  larger  slsed  orifice  holee  were  of  a  dlstinot 
advantage  because  of  the  greatly  reduced  responee  time  of  the  Instrument. 
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The  cro88~piece  for  both  probe8  is  at  a  sufficient  distance  downstream  from 
the  probe  head>  so  that  no  appreciable  flow  disturbance  occurs  at  the 
pickup, 

c,  Rotary-Seal  Pressure  Take-offs 

A  pressure  transfer  device  is  required  to  transfer  the 

pressure  readings  from  the  rotating  pitot-probe  to  stationary  manometers. 

12 

Hamrick,  Mizisin,  and  Michel  used  sealed  ball-bearing  units  for  the 
pressure-transfer  section.  Their  report  mentions  difficulties  incurred 
through  leakage  due  to  fai lures  of  the  seals.  As  a  preliminary  investiga¬ 
tion,  a  model  was  made  employing  various  commercial  brands  of  sealed  ball¬ 
bearings.  It  was  found  that  the  model  would  operate  perfectly  for  awhile, 
and  then  would  develop  intermittent  leaks.  In  view  of  the  uiureliable  per¬ 
formance  of  this  type  of  seal,  a  transfer  device  employing  contact-type 
rawhide  seals  was  developed.  This  seal  was  used  successfully  at  the  1500 
rpm  impeller  speed,  but  at  higher  speeds,  overheating  and  wear  became 
excessive  and  resulted  in  seal  failures.  Considerable  effort  was  eaqpended 
to  adequately  lubricate  aid  cool  the  seals,  and  at  the  same  time,  keep  the 
air  chambers  perfectly  free  from  the  lubricant.  Best  success  was  achieved 
using  air  cooling  and  small  amounts  of  molybdenum-disulfide  as  a  lubricant. 

In  an  effort  to  improve  the  pressure-transfer  device,  a 
viscous-drag  mercury  seal  was  designed.  A  pilot-model  with 
transparent  faces  was  made  to  determine  the  optimum  seal  gap,  mercury 
reservoir  volume,  and  range  of  operating  speeds.  Guided  by  the  tests 
conducted  with  this  model,  the  prototype  was  designed  and  is  shown  in 
Fig.  22  .  The  best  performance  was  obtained  by  using  a  smooth  disk  rotat- 
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Fig.  22  SECTION  OF  MERCURY  SEAL  FOR  PRESSURE  TRANSFER 
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Ing  in  the  fflereury  reservoir.  The  centrifugal  action  produced  by  disk 
friction  causes  the  mercury  to  rotate  and  fill  the  space  between  the  rotate 
ing  disk  and  the  stationary  housing.  A  disk  of  this  type  loeated  at  esoh 
side  of  the  pressure  take-off  opening  in  the  shaft  produces  a  sealed  air 
chamber  from  -hich  the  pressure  may  be  measured  with  a  Stationary  manometer. 

To  obtain  reliable  measurements  of  the  relative  flow  it  is 
necessary  that  the  sealing  be  perfect.  Errors  were  encountered  at  first 
because  of  contamination  of  the  mercury  caused  by  contact  with  base  metalsy 
0-rings,  or  minute  quantities  of  lubricants.  This  contamination  produces 
an  amalgamation  which  will  adhere  to  the  rotating  disk  and  Interfer  with 
the  cohesive  forces  between  the  d  Isk  and  the  mercury;  x^sulting  in  a  per-^ 
iodic  breakdown  of  the  sealing  surface  and  erratic  pressure  readings.  At 
speeds  greater  than  1000  rpm  the  seals  will  quickly  heat  up,  accelerating 
the  formation  of  mercuric  amalgams.  Air  cooling  was  used  at  the  1500  rpm 
speed,  but  was  Inadequate  for  higher  impeller  speeds.  Liquid  cooling  pre¬ 
sents  a  problem  because  of  the  possibility  of  very  small  amounts  of  the 
liquid  coming  Into  contact  with  the  mercury.  In  place  of  oil,  ethylene¬ 
glycol  was  used  as  the  coolant  because  it  has  less  serious  eontamlnation 
effects  on  the  mercury.  Contact  of  the  mercury  with  0-rings  tends  to 
promote  the  formation  of  amalgams  and  it  was  necessary  to  insert  a  tsflon 
back-up  ring  into  the  0-ring  groove  to  provide  a  barrier  between  the 
0-ring  and  the  mercury. 

3.  Overall  Performance 

The  pressures  at  the  inlet  filter  box,  and  Inside  the  discharge 
housing,  were  measured  to  determine  the  overall  energy  input  to  the  fluid. 
The  pressure  distribution  relative  to  the  blades  was  measured  at  the 
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impellar  outlet  by  means  of  static  pressure  taps  in  the  hub  surface  lead- 
ing  to  the  rotary  seals.  Static  pressure  taps  were  also  located  In  the 
walls  of  the  stationary  diffuser.  The  brake  torq^ue  and  power  were  measured 
with  the  dynamometer. 

C.  Test  Procedure 

The  velocity  distribution  relative  to  the  rotating  iii5)eller  is 
measured  at  the  blade  outlet  with  the  spherical  pitot-probe  mounted  at  the 
six  circumferential  positions  shown  in  Fig.  23.  An  eight-point  traverse 
(Tig.  2h  j  is  Bude  at  each  of  these  six  circumferential  positions.  The  re^ 
suiting  U3  positions  cover  the  discharge  area  of  1^.7  square  inches  be¬ 
tween  two  adjacent  blades.  Two  separate  runs  were  required  at  each  test 
point  in  order  to  check  for  leaks.  This  was  accomplished  by  interchanging 
the  tubes  leading  from  the  pitot-probe  holes  to  the  rotary  seals,  and 
checking  the  reproducibility  of  the  manometer  readings,  A  total  set  of 
readings  were  taken  for  a  single  speed  and  a  fixed  flow  rate.  The  flow 
rate  was  measured  at  the  standard  inlet  nozzle  and  adjusted  by  means  of 
the  louvers  at  the  outlet  housing.  The  pressures  in  the  inlet  filter 
chamber,  the  inlet  draft-tube,  the  diffuser  walls,  and  the  drum  housing 
at  the  outlet  were  measured.  After  the  internal  measurements  of  the  flow 
distribution  were  completed,  the  rotary  seal  assembly  was  removed  to 
eliminate  the  torque  of  this  unit,  and  the  overall  performance  of  the 
iaqseller  was  measured  with  the  dynamometer  drive. 

P.  Presentation  of  Data 

1.  Internal  Flow  Measurements 

The  results  of  the  detailed  measurements  of  the  flow  distribu¬ 
tion  at  the  impeller  outlet  are  presented  in  Table  V.  The  first  sections 
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of  the  table,  (a)  through  (g),  conQirise  the  data  taken  with  the  spherical 
probe  with  the  larger  orifice  holes  (0.016  inch  diameter),  and  using  the 
mercury  seals  in  the  pressure  transfer  device*  The  latter  sections  of  the 
table,  (h)  through  (k),  present  the  data  taken  with  the  spherical  probe 
with  the  smaller  orifice  holes  (O.Oll  inch  diameter)  and  Slower  response 
time,  and  using  the  contact- type  rawhide  seals  in  the  pressure-transfer 
device.  The  discharge  rate  is  considerably  lower  for  the  runs  made  with 
the  rawhide  seals  because  a  small  high-resistance  filter  was  used  at  the 
inlet  to  the  draft  tubej  also,  for  these  runs,  a  straight  conical  inlet 
to  the  draft-tube  without  straightening  vanes  was  used  in  place  of  ^ 
bell-mouth  nozzle. 

Graphs  of  the  data  showing  the  distribution  of  the  pertinent 
flow  variables  are  given  . in  Figs.  25  through  30  .  The  distribution  of  the 
meridional  coiqjonent  of  the  relative  velocity  across  the  channel  from  hub 
to  shroud  is  shown  in  Fig.  25.  For  purposes  of  comparison,  the  theoretical 
velocity  distribution  for  the  meridional  flow  pattern,  undisturbed  by  the 
blades,  is  shown  for  each  of  the  flow  rates.  The  distribution  of  the 
meridional  coa$)onent  of  the  relative  velocity  from  blade  to  blade  is  Shown 
in  Fig.  26 ,  and  the  direction  of  the  relative  velocity  in  terms  of  d  and  p 
is  shown  in  Figs.  27  and  28.  Figures  29  and  30  Show  the  distribution  of 
the  whirl-moment  and  specific  energy.  In  Figs.  26  throu^  30,  the  plots 
for  the  data  taken  with  the  rawhide  Seals  are  shewn  only  for  the  mean 
meridional  streamlines. 

2.  Efficiency,  Dischargej  and  Head  Characteristic 

The  primary  objective  of  the  ej^rimental  investigation  was 
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FOUR-BIADED  IRPELLER  DATA  K>im-SUDED  1MPEU£R  DATA 
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to  obtain  neasurements  of  the  internal  flow  distributions.  The  test 
stand  was  designed  with  this  in  mindf  and  the  outlet  diffuser  is  eon> 
structed  so  that  disturbances  of  the  relative  flow  will  be  minimized. 

Tor  this  reason  a  collector  ring  is  not  included,  and  the  conventions^ 
compressor  efficiency  which  Includes  the  diffuser  performance  cannot  be 
measured.  However,  as  an  indication  of  the  general  performance  of  the 
impeller  design,  the  shaft  power  was  measured  with  the  dynamometer,  and  the 
corresponding  static  pressure  was  recorded  for  a  range  of  flow  rates  at 
3500  rpm.  An  "overall  efficiency"  was  computed  by  considering  the  energy 
input  to  the  gas  to  be  the  sum  of  the  static  pressure  head  in  the  diffuser^ 
plus  an  average  velocity  head  based  upon  the  flow  rate  and  the  diffuser 
area.  This  overall  efficiency  is  shown  in  Fig.  31 .  It  must  be  noted  that 
this  efficiency  includes  the  draft-tube  inlet  losses;  there  is  also  prob¬ 
ably  a  considerable  windage  loss  caused  by  the  stainless  steel  tubing. 
Supporting  ring;  and  pitot*positloning*di8k  on  the  back  of  the  impeller. 

The  same  energy  input  considered  for  the  determination  of  the  overall  ef« 
ficiency  was  also  used  for  the  head  characteristic.  The  head  input  is 
plotted  as  a  dimensionless  ratio  based  upon  the  impeller  tip  speed.  The 
impeller  efficiency  is  determined  for  the  2200  rpm  speed  from  the  data 
presente<l  in  Table  V.  The  useful  energy  H  of  the  gas  is  considered  to  be 
the  sum  of  the  average  static  pressure  head  and  the  velocity  head,  and  the 
impeller  work  is  com/g. 

The  total  flow  entering  the  draft  tube  is  measured  with  the 
standard  A.S.H.E.  long-radius  nozzle.  The  A.S.M.E.  Fluid  Keber  Report  was 
used  to  determine  the  variation  of  the  discharge  eoefficieht  with  the 
Reynolds  number.  The  discharge  was  also  confuted  at  the  impeller  outlet 
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by  a  numerical  integration  of  the  meridional  velocity  distribution,  A 
comparison  of  the  rates  of  flow,  as  measured  with  the  nozzle  and  computed 
from  the  velocity  survey,  is  shown  in  Table  VI, 

IV.  GORRELATION  GF  THEORY  AND  EXPERIMENTAL  RESULTS 


•  Msturbance  of  Meridional  Streamline  Pattern 


The  impeller  design  is  based  upon  the  assumption  that  at  the  design 
Condition  of  Q/(oD^  ■  0,060,  the  meridional  velocity  distribution  is  un¬ 
disturbed  by  the  presence  of  the  blades.  The  distribution  of  the  theoret¬ 
ical  velocity  shows  a  uniform  increase  from  the  hub  to  the  shroud.  The 
Velocity  at  the  hub  is  13  per  cent  below  the  average,  and  at  the  Shroud  the 
Velocity  is  16  per  cent  above  the  average.  Ebcamination  of  Pig,  25  shows 
that  this  general  trend  for  the  meridional  velocity  to  increase  uniformly 
from  hub  to  shroud  occurs  for  all  of  the  test  rxins,  and  for  values  of  the 
discharge  parameter  Q/wD^  ranging  from  0,023  to  0,087,  However,  while  the 
shape  of  the  curve  does  not  significantly  change,  it  is  displaced  in  the 
peripheral  direction  about  ±  10  per  cent  from  the  average  (theoretical 
distribution),  with  the  exception  of  the  c\urve  immediately  behind  the  trail¬ 
ing  edge  of  the  hlade  tip  (0  >■  82°  20* ) .  The  meridional  velocity  dietributtaa 
just  behind  the  trailing  edge  is  displaced  about  20  per  cent  above  the  aver¬ 
age,  and  the  reduction  in  the  velocity  profile  due  to  boxuadary  layer  build¬ 
up  starts  at  a  considerable  distance  from  the  shroud*  This  effect  becomes 
more  pronounced  at  the  lover  discharge  rates,  indicating  a  significant 
boundary  layer  buildup  . at  the  shroud  surface  near  the  tarailing  edge.  The 
large  pressure  gradient  across  the  channel  from  hub  to  shroud,  occurring 
behind  the  trailing  edge,  may  indicate  lui  «g>preciable  secpndury  flow  in 
this  region,  (Figs.  3^  and  33)* 


32  STATIC  FSESSOHE  VAKXATIOR  FROM  HUB  TO  SHROUD  AT  MPMJJtR  OUflAT 


0  30*  60*  90^ 

Angular  Displacenent  From  Leading  Edge  of  Blade >  Q 


Fig.  33  STATIC  PRESSURE  VARIATION  PROM  BLADE  TO  BLADE  ON 
STREAM  SURFACE  AT  IMPELIIR  0UT1£T. 
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B.  Similarity  Considerations 

It  ia  generally  accepted  that  similitude  relations  are  satisfied 
for  turbomachines  operating  under  the  conditions  of  this  investigation 
(Vislicenus^^,  Stanlts^)  in  which  Hach  number  and  Reynolds  number  effects 
are  not  involved.  This  means  that  the  velocity  diagrams  at  all  points  of 
the  flow  field  will  remain  similar  over  a  range  of  speeds,  provided  the 
discharge  parameter  is  maintained  constant.  Comparison  of  the  di¬ 

mensionless  velocity  distributions  of  Figs.  25  and  26,  and  of  the  velocity 
direct! on&  shown  in  Figs.  27  and  28,  shows  that  similitude  relations  are 
effectively  maintained,  and  tends  to  verify  the  reliability  of  the  measure¬ 
ments.  Because  of  this  similarity,  it  is  sufficient  to  concentrate  on  the 
data  for  one  particular  speed,  2200  rpm,  to  study  the  variations  occuring 
at  the  design  conditions.  Complete  data  at  this  speed  are  presented  for 
a  range  of  values  of  of  0.023  ,  0.060,  and  0.085>  which  bracket  the 

design  value  of  0.060. 

C.  Variation  in  the  Energy  and  Vfliirl«^Moment  Values 

Icccardlng  to  the  infinite  blade  theory,  the  energy  input,  and  the 
increase  in  the  whirl’^moment,  should  be  uniformly  distributed  at  the  blade 
outlet.  The  variation  in  these  values  is  shown  in  Figs.  29  and  30.  In¬ 
spection  of  the  (c)  and  (d)  parts  of  these  figures  for  the  2200  rpm  speed, 
just  above  and  below  the  design  condition,  reveals  a  small  variation  of  * 
10  per  cent  (with  the  exception  of  the  trailing  edge  of  the  blade)  across 
the  channel  from  hub  to  shroud.  However,  there  is  a  significant  variation 

from  blade  to  blade.  Ibis  same  type  of  variation  has  been  observed  by 
12 

Hamrick,  et  al.  .  The  maximum  energy  input  occurs  at  the  driving  face 
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of  the  bLadef  supporting  the  hypotheels  made  In  Eeferenee  12,  that 
"eeeondary  fLows  In  the  boundarj  layer  tend  to  ehift  the  low«ener|Qr  «ir 
toward  the  blade  suction  surface"* 

The  ralue  of  1  should  remain  constant  for  isentropic  flow,  and  is 
a  measure  of  the  "loss",  which  is  defined  as  the  difference  between  the 
"nseful"  energy  H  and  the  energy  input  to  the  gas  given  by  dm/g.  The 
iaqpeller  efficiency  is  defined  as  the  ratio  of  these  two  values.  The  vaf* 
iatlon  in  the  impeller  efficient  is  shown  in  Fig.  3it-.  These  efficiency 
values  include  the  eneror  losses  in  tiie  draft  tube.  The  lowest  efficiency 
occurs  near  the  shroud  surface,  across  the  entire  passage  from  blade  to 
blade;  the  losses  are  particularly  high  in  the  corner  fcrmed  by  the  trail> 
ing  face  of  the  blade  and  the  shroud  surface. 

D.  Comparison  with  the  Blade  to  Blade  Solution 

For  the  blade  to  blade  solution,  it  is  assumed  that  the  flow  pro^ 
ceeds  along  surfaces  of  revolution  generated  by  meridional  streamlines* 
Variation  of  the  flow  angle  d  along  the  periphery  is  a  measure  of  the  ex¬ 
tent  bhat  the  flow  "twists"  away  from  this  surface  of  revolution.  Exam¬ 
ination  of  Fig,  29d  (2200  rpm,  ■  0.060)  shows  that  the  angle  8 

varies  only  about  3  degrees  on  any  given  geometric  surface  of  revolution. 
This  indicates  that  the  fluid  particles  remain  essentially  on  surfaces  of 
revolution  for  the  main  flow  region  between  the  blades.  The  uniformity 
of  the  flow  from  hub  to  shroud  indicates  that  a  satisfactory  picture  of 
the  flow  across  the  tfitire  channel  will  be  derived  from  values  co)q}uted 
for  the  mean  stream  surface. 

Comparison  of  the  theoretical  meridional  velocity  distributicus 
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with  th«  MtsTirad  retulti  1#  Bho«n  in  Fig.  3^*  Xt  ma^  observed  that 
the  general  trend  in  the  measured  distribution  is  similar  to  the  theoreti¬ 
cal  for  all  values  of  the  discharge  parameter.  The  measured  values  do  not 
rise  as  steeply  as  the  theoretical  near  the  trailing  face  of  the  blade. 

This  is  most  likely  caused  by  unloading  of  the  pressure  forces  on  the 
blade  near  the  blade  tip,  because  the  pressure  difference  across  the  blade 
must  become  aero  at  the  tip.  This  hypothesis  is  supported  by  the  distribu¬ 
tion  of  static  pressure  from  blade  to  blade >  as  shown  in  Fig«  33.  The 
static  pressure  is  observed  to  be  a  minimum  a  short  distance  from  the  suc¬ 
tion  side  of  the  blade}  and  then  increases  upon  approaching  the  blade. 

This  indicates  the  presence  of  a  secondary  flou}  tending  to  produce  a 
velocity  distribution  at  the  blade  outlet,  which  is  more  uniform  than  that 
computed  on  the  basis  of  potential  theory. 

1.  The  Iiig)eller  Slip  Factor 

The  slip  factor  /4  is  an  is^ortant  parameter  for  the  determination 
of  the  ii^ller  performance,  and  is  defined  as  the  ratio  of  the  average 
value  of  the  tangential  co^)onent  of  the  absolute  velocity  divided  by  the 
theoretical  value  that  would  be  obtained  if  the  fluid  were  perfectly 
guided  by  the  vanes.  Stodola^^  gives  an  equation  for  straight  blades, 

fi-  1  -  0.500  AO  , 

and,  for  a  blade  spacing  of  1.502  radians,  the  slip  factor  Should  be 
0.2U9.  Stanits^  is  not  in  agreement  with  the  Stodola  correction;  he 
has  computed  the  slip  factor  by  relaxation  methods  for  a  series  of  eentri- 
fugal  impellers  having  straii^t  blades,  and  also  for  two  impellers  having 
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0  30*  6G*  90* 

Angular  Dlaplacement  From  Leading  Edge  of  Blade >  0 


Fig.  33  COMPARISON  OF  KEASNREO  AND  THEORETICAL  MERIDIONAL  VELOCITY 
DISTRIBUTION  ON  MEAN  STREAM  SURFACE  AT  THE  SfFELLER  OUTUT. 
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logaflthmle-splral  blades  of  backward  curvature,  and  found 
^  -  1  -  0.315a©  . 

This  equation  leads  to  a  slip  factor  of  0,52?  for  the  foir-blade  impeller 
of  this  investigation.  A  comparison  of  the  measured  slip  factors  is 
shown  in  Fig.  36,  This  figiure  indicates  that  the  slip  factor  depends 
upon  the  speed  and  discharge  in  an  actual  impeller,  and  shows  a  var¬ 
iation  of  nearly  ±  12  per  cent  about  a  median  value  of  0.571,  leading  to 
an  equation  for  the  slip  factor, 

^  ■  1  -  0.30A6  , 

which  is  in  close  agreement  with  the  theoretical  relationship  of  Stanitz. 

F,  Overall  Performance 

The  large  deviation  between  the  impeller  efficiency  and  the  over¬ 
all  efficiency  shown  in  Fig.  31,  particularly  at  low  flow  rates,  indicates 

* 

the  existence  of  eddy  motion  within  the  impeller.  The  existence  of  an  eddy 
attached  to  the  driving  face  of  the  impeller  at  low  flow  rates  is  theoret¬ 
ically  shown  by  Stanitz^^,  Reverse  flow  on  the  driving  face  is  also  in- 

17 

dicated  for  flow  solutions  at  low  weight  flows  by  Prian  and  Michel  but 
direct  experimental  evidence  of  the  eddy  was  not  found.  A  possible  ex* 
planation,  offered  in  reference  (12),  Is  that  the  flow  is  unstable,  and 
that  instantaneous  eddies  will  form  and  be  washed  awey, somewhat  in  the 
manner  of  rotating  stall  observed  in  axial-flow  compressors.  Detection 
of  this  phenomena  is  not  possible  with  the  spherical  pitot-probe,  as 
this  instrument  will  only  give  a  time  average  of  the  total  pressure,  and 

*  See  reference  (12) 


WADC  ra  55-158 


79 


this  tliM-*Terage  Is  impossible  to  interpret  in  areas  of  rapidly  changing 
flow  direction*  This  instability  mny  explain  the  discrepancy  in  the  total 
floH  rate  as  meastu>ed  at  the  inlet  noszle,  and  computed  by  numeric^  in¬ 
tegration  at  the  i^)eller  outlet,  shoirn  in  Table  VI,  at  the  low  rate  of 
flow.  The  integrated  value  is  nearly  27  per  cent  lower  than  ^  value 
measured  at  the  nozzle, 

V,  SUIUARY  AND  GQWGLUSIOMS 

A  four-blade, mixed- flow  impeller  has  been  designed  on  the  basis  of 
axisymmetric  flow  conditions  derived  from  the  assueption  of  an  infinite 
number  of  blades.  The  flow  distribution  at  the  outlet  of  this  impeller 
has  bsen  measured  with  a  spherical  pitot-probe  mounted  on  the  impeller,  and 
the  results  of  these  measurements  are  presented  in  graphical  form.  Compar¬ 
isons  between  the  theoretical  flow  distributions  and  tha  measured  distribu¬ 
tions  have  been  made,  and  the  following  conclusions  are  drawn. 

1.  The  shape  of  the  meridional  velocity  distribution  remains  essen¬ 
tially  unchanged  from  the  theoretically  predicted  shape  based  upon  the  as- 
sun^tlon  Of  axial  symmetry,  in  spite  of  the  small  blade  number  (four)  of 
the  tBq3eller.  However,  in  tha  blade  to  blade  direction,  the  meridional 
distribution  is  displaced  as  much  as  20^1  above  and  IO56  below  the  aver¬ 
age,  This  trend  is  observed  for  a  range  of  the  discharge  parameter 
from  0.023  to  O.Od^. 

2*  The  blade  surface  Is  considered  to  be  a  bound -vortex  sheet  with  a 
zero  ring-component  of  vorticlty.  This  means  the  blades  do  not  induce 
velocity  components  which  disturb  the  meridional  velocity  distribution  at 
the  design  condition.  The  observed  conformity  of  the  meridional  velocity 
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dlatrlbutlon  with  the  theoretical  distribution  for  potential  through-flow 
in  the  annular  space  formed  by  the  hub  and  shroud  0  urfaces^  tends  to  verify 
the  vortex  theory  of  the  flow. 

3*  the  fluid  particles  remain  essentially  on  surfaces  of  revolution 
generated  by  the  meridional  streamlines.  Variations  in  the  cone  angle  S 
of  the  relative  velocity  are  of  the  order  of  only  2  and  3  degrees  in  the 
peripheral  direction.  This  substantiates  the  underlying  assumption  of 
the  quad i^ three-dimensional  relaxation  solutions  of  the  flow. 

U*  the  observed  blade-to-blade  velocity  distribution  is  more  uniform 
than  the  distribution  computed  by  relaxation  methods.  This  may  be  due  to 
the  large  blade  spacingi  and  the  necessary  blade  unloading  near  the  tip 
which  tends  to  increase  the  pressure  at  the  suction  (trailing)  face  of  the 
blade}  and  decrease  the  pressure  on  the  driving  (leading)  face  of  the 
blade.  This  unloading  tends  to  smooth  out  the  abrupt  velocity  change 
across  the  blade}  as  is  indicated  by  the  relaxation  solution. 

Low  energy  air  tends  to  shift  towards  the  suction  side  of  the  bladeS} 

12  17 

which  coincide  with  observations  made  by  other  investigators  ’  . 

This  shift  is  apparently  caused  by  secondary  flows  in  the  boundary  layer} 
and  may  account  for  the  discrepancy  in  the  relaxation  solution  near  the 
blade  surface. 

6.  The  Impeller  slip  factor  varies  nearly  12  per  cent  about  a  median 
value  of  0.57}  over  a  range  of  speeds  (1500  to  3000  rpm)  and  discharge 
rates  (qA)D^  •  0.023  to  0.085)  •  The  measiured  slip  factor  compares  favor¬ 
ably  with  a  theoretical  value  of  0.53  given  by  StanltS}  but  is  In  disagree¬ 
ment  with  the  value  of  0.25  eoiqputed  f  jram  ^e  Stodola  equation  for  centrl- 
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fugal  impellers  having  straight  blades. 

Regarding  the  experimental  technique  applied  during  the  investigation, 
it  can  be  concluded  that  consistent  and  reproducible  surveys  of  the  internal 
flow  distribution  can  be  obtained  with  a  directional  pitot-probe  mounted 
on  the  impeller  and  used  in  conjunction  with  a  mercury-seal  type  of 
pressure-transfer  device.  The  readings  become  unreliable  if  flow  instabil¬ 
ities  resulting  from  separation  and  eddy  motion  are  present. 
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APPENDIX  A 
SYMBOLS 

I.  COQRDINATB  SYSTEMS 

(r,0,z;  Cylindrical  coordinates  rotating  with  the  impeller 

(RiOyZ)  Dimensionless  cylindrical  coordinates,  where 

R  ■  r/D,  Z  ■  z/D 

(8,n,0)  Meridional  coordinates  indicating  the  arc-length  along 
a  meridional  streamline,  the  normal  to  the  meridional 
streamline  (positive  towards  the  center  of  curvature), 
and  the  angular  location  around  the  impeller  axis  of 
rotation. 

(S,N,Q}  Dimensionless  meridional  coordinates,  where 

S  -  s/D,  N  -  n/D 

Orthogonal  coordinates  on  the  mean  stream  surface  de¬ 
fined  by  equation  39. 

n.  VECTORS 

Vectors  are  indicated  by  a  bar  over  the  letter.  The  vector  product  is 

indicated  by  a  cross  x  ,  and  the  sealar  product  by  a  dot  •  . 

III.  VELOCITIES 

“c  Absolute  velocity  of  a  fluid  particle,  c  »  u  ♦  w 

u  Peripheral  velocity  of  a  point  of  the  rotating  impeller, 

u  •  rw 

w  Velocity  of  a  fluid  particle  relative  to  the  impeller. 

(*)  Angular  speed  of  the  impeller. 

C,W  Dimensionless  velocities,  C  ■ 

IV.  OPERATOB3 

e3q>  Exponential,  (exp  x  ■  e*) 

7  Gradient  vector  operator 

7^  Laplacian  operator,  (7*7) 
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V.  INDimS 


Differentiation  following  the  motion  of  the  fluids 


Small  Incremental  variation 


Increment 


r«t|a  Represent  components  in  the  radial,  tangential  and  axial 

directions  respectively 

8,n  Represent  components  in  the  meridional  coordinate  $ystem 

m  Component  in  the  meridional  plane 

rel  Relative  to  the  blade  cascade 

(other  indices  used  are  explained  in  the  text,  or  are 
considered  an  part  of  the  definition  of  the  symbol.) 

VI,  SPECIAL  NOTATIONS  FOR  THE  SiH^JilCAL  PITOlWPHQBE  (See  Fig,  Bl) 

(x',y*»*')  Rectangular  coordinates  fixed  to  spherical  pitot-probe 

with  x’-axis  along  the  sphere  axis,  and  y'-axis  parallel 
to  the  cross-piece. 

Unit  vectors  in  the  (x*,y',z')  system 


Unit  vectors  in  the  (x*,y',z')  system 

e  ,e.,e  Unit  vectors  in  the  (r,0,z)  system 

n  V  s 

a  Conical  angle  formed  by  the  four  synmetrical  pressure 

holes  and  the  probe-axis  (x'-axis). 

P  Angle  between  the  probe-axis  (x'-axis)  end  the  tangential 

^  direction  (t-dtrection) 

y  Angle  between  the  relative  velocity  vector  and  any  pres- 

”  sure  hole  (n)  on  the  spherical  probe-head. 

5'  Dihedral  angle  between  the  plane  containing  the  incident 

relative  velocity  and  the  meridional  plane  of  the 
spherical  probe-head, 

0  Conical  angle  formed  by  the  relative  velocity  and  the 

probe-axis  (x'-axis) 

cos(x',r)  The  direction  cosine  between  the  x'  and  r-axis,  etc. 
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Vli.  OTHERS 


A  Acceleration  vector  of  a  fluid  particle 

6  Dimensionless  height  of  axial-symmetry  flow  channel 

formed  by  stream  Surfaces  adjacent  to  the  surface  of 
revolution  generated  by  the  mean  meridional  streamline^ 

B  - 


c 


mo 


Absolute  meridional  velocity  component  at  the  hub  surface 
(n-o) 


D  Maximum  diameter  at  blade  outlet-edge 

F(r,z)  Fvinction  defining  the  blade  surface 
T(0»8(0  Functions  defined  by  equations  (li5)  and  (li6) 


8 

H 

h 

s 

h 

I 

k 

m 

n 

An 

P 

Q 

<3n 


Acceleration  of  gravity 

Total  head  input  to  fluid  (specific  energy),  H 
Static  pressure  head 


2g  i  J  f 


Relaxation  grid  spacing 
u 

i 


Fictive  distributed  blade  force  per  unit  volume  (for 
infinite  number  of  blades) 

Whirl-moment,  m  •  rc^ 


Normal  vector  to  the  blade  surface 

Spacing  between  hub  and  shroud  surfaces  along  orthogonal 
in  the  meridional  plane 

Pressure 


Volume  rate  of  flow  through  impeller 


Volume  rate  of  flow  between  the  hub  surface  and  the 
Stream  surface  of  revolution  generated  by  meridional 
streamline  (n) 

Residual  (error  in  finite-difference  equation  due  to 
incorrectly  estimated  function  values.) 
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Angle  between  relative  fluid  velocity  and  negative  blade 
velocity,  between  w  and  (“u) 

Circulation 

Cone  angle  formed  by  direction  of  meridional  component  of 
relative  fluid  velocity  with  impeller  axis  (between  W- 
and  z)  ”* 

Impeller  efficiency,  ■  glV^m 

Absolute  vorticity 
Absolute  bound  vorticity 

Absolute  free  vorticity 
Impeller  slip  factor ,^4  •  (C^) 

average  theoretical 

Fluid  mass  per  unit  volume 

Stream-^f unction  defined  by  equation  (38) 

Change  in  \j/  at  grid  points 

Dimensionless  stream^function  indicating  the  quantity  of 
flow  passing  between  the  ’'I'-streamline  and  the  leading 
face  of  the  blade,  expressed  as  a  percentage  of  the  flow 
between  two  blades^  IqO  JiL 


t 


APPENDIX  B 

THE  THREE  DIMENSIONAL  SPHERICAL  PITOT  PROBE 

1,  Theoretical  Calibration  Curves 

The  theoretical  relationship  between  the  relative  velocity  of  the 
free*stream  at  a  point,  and  the  pressure  readings  given  try  the  five  holes 
of  the  spherical  probe-head  of  Fig.  (Bl)  if  it  is  inserted  into  the  stream 
at  this  given  point,  will  be  derived.  As  usual,  the  fluid  will  be  assumed 
to  be  nonvisGoue,  incompressible^  and  steady  relative  to  the  probe,  so 
that  the  Bernoulli  Equation  will  apply  between  a  point  in  the  free  stream 
and  hole  (n)  on  the  ^}here, 

«  ..  2 


P  •  Pp  ♦  /  n  •  0,  1,  2,  3»  U. 

Solving  this  equation  for  the  pressure  at  hole  (n)  yields, 

2 


(Bl) 


Pn  -  P  ♦PT 


•  P  *  ‘'n^T 

The  pressure  recovery  factor,  k 


a  function  of  the 


rw 

n  "  V  w . 

Reynolds  nuntrer  and  the  orientation,  given  by  the  angle  of  the  hole 
with  respect  to  the  stream.  However,  this  factor  becomes  independent  of 
the  Reynolds  n\imber  fcr  sufficiently  high  velocities,  and  is  only  a 
fxinction  of  which  can  be  expressed  in  terms  of  the  two  angles  0  and  S'. 
By  combining  the  equations  for  the  pressure  at  each  of  the  five  holes, 
the  following  ratios  may  be  obtained: 
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•f 


P0-P3 


(B3) 

(BU) 


Each  of  these  ratios  may  be  determined  experimentally  by  measuring  the 
pressure  differences  occurring  between  the  corresponding  holes  on  the 
pitot  sphere  for  a  known  orientation  of  the  probe  with  respect  to  ^e 
free-Stream  velocity  vector.  The  two  functions  of  equations  (B3)  and  (Bl*) 
implicitly  give  0  and  6*,  which  determine  the  direction  of  the  velocity, 
in  terms  of  the  ratios,  (k2-k^)/(k^-k^)  and  (k^*i-k^)/(k^**k^) . 

These  ratios  have  finite  values  if  5*  is  in  the  range  0<6'  <180* 
(Fig.  Bl).  But  in  the  range  -180*<6'  <0,  the  ratios  become  infinite  when 
the  vector  makes  equal  angles  with  holes  0  and  3»  For  this  reason,  and 
also  for  the  sake  of  symmetry,  the  graph  of  the  two  ratios  defining  flow 
direction  and  velocity  is  divided  into  two  halves,  0<6'<180*,  and 
-180*<5’<0,  For  the  first  half,  the  ratios  given  by  equations  (B3)  and 
(Bit)  are  used.  For  the  second  half,  the  ratios  used  arei 


k,-k  p,-p 


‘o’"!  Po’Pl 


(B7) 


(B8) 


The  theoretical  determination  of  the  above  ratios  of  the  pressure 
recovery  factors  may  be  based  upon  potential  flow  around  a  sphere.  For 
this  case,  the  velocity  at  any  point  on  the  sphere  surface,  corresponding 
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to  a  radial  line  forming  an  angle  with  the  wadirectiony  is  given  by 
the  equation 


Wn  -  (|  sin  yj  w  . 

This  means  that  the  pressure  recovery  factor  at  hole  (n)  for  potential 
flow  is  given  ly 


k  •  1  -  ^  ain^  y  , 
n  a  '  n' 


and  the  ratio  of  the  pressure  recovery  factors  is 


k^-k. 


2  2 
sin  )c-sin  y^ 

b"*‘3  sin^y^-sin^yQ 


2  2 
ooa  yj^-cos  y^ 

“  — 2" - T- 

cos  J^’^cos  y^ 

The  ratios  of  the  other  pressure  recovery  factors  are  given  by  similar 
expressions • 

The  angle^j^  is  related  to  0  and  6'  by  the  equations  of  spherical 
trigonometry, 

cos  y^  ■  cos  0  , 

cos  y^  -  cos  0  COSO  ♦  sinOsinS'sin  0  , 

cos  y2  •  cos  0  COSO  sinocosS'sin  0  , 

cos  y^  ■  cos  0  COSO  ^  sinosinS'sin  0  , 

cos  y^  ■  COS  0  cosO  •-  sinacos6'sin  0  • 


(filO) 
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Kaking  use  of  these  relations,  the  functions  representing  the  theoretical 
ratios  of  the  pressure  recovery  factors  are 


2  slh  2a  sindVsin  2(jl( 

2  3in^a(cos^0-  sin^S'  3in^0)-*.  sin  2a  sinS'sin  20 
2  sin  2a  cos  S*  sin  20 

2  sin^  (cos^0-  sin^S*  sin^0)>  sin  2^  sin  S' sin  20 


(B12) 


These  ratios  are  plotted  in  Fig,  B2.  Because  of  perfect  symmetry,  only 
the  first  quadrant,  in  which  S'  ranges  from  zero  to  90  degrees  is  shown. 
The  above  functions  are  used  to  determine  the  direction  o  f  the  vel¬ 
ocity,  The  magnitude  of  the  velocity  may  be  determined  by  applying  the 
general  pressure  equation  to  the  center  hole,  and  to  any  one  of  the  other 
four  holes,  say  hole  (3),  which  will  give  the  equation 

^2 

VP3  “^T  ^V^3^* 

From  this,  the  velocity  is 


Equation  (B8)  and  (BlO)  may  be  used  to  determine  the  theoretical  relation 
of  (k^-?k^)  to  0  andS'  , 


VS  “  f(cos^y^-coaV3) 

■  ^  jsin^a  (cos^0^  ain^S’  sin^0)+  ^  sin  ZQ  sin  20  sin  6'  (BlU) 

After  the  magnitude  and  direction  of  the  relative  velocity  has  been 
determined,  the  static  pressure  of  the  free  stream  is  given  directly  by 
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^^*?^'*'’*^g^'^W^S^!W*^'^*l•'y■!^f*»^'^*^>l^Vt*^^^  -V  • '  *•  •■  <piiir.w^<u»«i  up'iyfi’gWBiWMl'l 


applying  equation  (B2)  to  the  center  hole« 

w2 

P  "  P©  ”  *  (®i5) 

The  theoretical  value  of  is  given  by  equation  (B8),  where  is 
Identical  to  the  cone  angle  0, 

ko  *  1  *  I  ain^  0  .  (Bl6) 

2^  Experimental  Calibration  Curves 

The  spherical  probe  was  calibrated  in  .the  test  stand  shown 
schematically  in  Fig.  B3.  The  calibration  fixture  consists  of  two  concen¬ 
tric  rings  in  a  glmble  mounting.  The  bottom  ring  is  pivqted  along  its 
diameter  in  the  table-mounting,  with  the  pivot-axis  perpendicular  to  the 
nozzle-axis.  The  tilting  angle  is  measured  with  a  protractor  on  the 
pivot-axis,  and  is  equal  to  the  cone  angle  0  between  the  velocity  vector 
of  the  flow  and  the  shaft-axis  of  the  probe.  The  top  ring  rotates  over 
the  bottom  ring,  which  has  protracted  angular  divisions,  to  measure  the 
dihedral  angle  5'.  The  probe  is  mounted  on  the  top  ring  with  its  axis 
coinciding  with  the  ring-axis,  and  with  the  center  of  the  spherical  probe 
on  the  pivot-axis  of  the  bottom  ring.  By  this  arrangement,  the  center  of 
the  sphere  remains  on  the  nozzle-axis  and  fixed  in  position,  for  aiy  change 
of  cone  angle  aid  dihedral  angle.  Thus,  as  these  angles  are  changed,  the 
error  caused  by  non-uniformity  of  the  velocity  distribution  of  the  air 
emerging  from  the  nozzle  is  minimized. 

Two  spherical  probes  were  designed  and  calibrated.  The  construc¬ 
tional  details  of  the  first  probe  that  was  used  are  shown  in  Fig.  iUa. 
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Fig.  B3  SPHERICAL  PITOT-FROEE  CALIBRATION  FIXTURE. 
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the  inside  diameter  of  the  hypoderndc-needle  tubing  is  0.011  inch,  the 
high  resistance  caused  hy  this  rather  long  length  of  small  tubing  resulted 
in  a  slow  response  time  for  the  instrument.  A  second  probe  was  designed 
having  a  shorter  length  of  small-diameter  tubing,  and  the  Inside  diameter 
was  increased  from  0.011  inch  to  0.016  inch,  the  details  of  this  second 
probe  are  shown  in  Fig.  BUb.  the  original  probe  could  be  locked  in  only 
one  position,  which  was  not  adequate  to  cover  the  wide  variations  encountered 
in  the  exit  flow  angles.  Qulded  by  the  measurements  obtained  with  the 
original  probe,  two  flats  were  ground  on  the  crosspiece  of  the  second  probe 
so  that  it  could  be  locked  in  the  two  most  favorable  positions. 

The  experimental  calibration  curves  for  the  original  pitot-probe 
are  shown  in  Flgs.B^a  and  B6a,and  for  the  second  modified  probe  In  Figs. 
B5bandfl6b.  The  general  trend  of  these  e:cperimental  curves  conforms  to  the 
theoretical  curves,  more  closely  at  small  angles  of  0,  and  for  increasing 
angles,  the  e3q)erimental  curves  become  markedly  different  from  the  theo¬ 
retical,  By  coo^aring  the  calibration  of  the  two  different  probes,  one 
may  derive  a  measure  of  the  individual  discrepancies  that  can  be  expected 
in  probes  that  have  the  same  specified  dimensions.  On  the  basis  of  this 
comparison,  if  flow  avgles  are  to  be  measured  within  an  accuracy  of  the 
order  of  one  or  two  degrees,  it  is  necessazy  to  calibrate  each  instriunent 
separately. 

3.  Relation  of  the  Probe  to  the  Inpe Her  Geometry 

The  calibration  curves  of  the  spherical  pitot«*probe  will  give  the 
magnitude  and  direction  of  the  fluid  velocity  relative  to  the  probe.  It 
is  necessary  to  express  this  measured  velocity  in  relation  to  the  blade 
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(b)  Orifice  Diameter  O.OI6  inch 


Fig.  B6  CALIBBATIGN  OF  SEHEEICAL  PITOT  TOR  DETERMmTION  OF  STXTIE  PHESSUBB, 
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cascade.  To  do  this  requires  a  transformation  from  the  coordinates 
(x',y',z')  fixed  to  the  pitot-probe,  to  the  impeller  coordinates  (r,Q,z). 
Referring  to  Fig.  B1  in  conjunction  with  Fig.  2k,  the  following  direction 

cosines  are  immediately  evident: 

cos(y*,z)  ■  -  cos  21*12'  ■  -  0.9323* 
cos  (y',r)  -  sin  21*12'  -  O.36I6, 
cos(y',t)  -  cos  90*  *  0, 
cos(x',t)  *  cosj3p, 
cos(z',t)  ■  sin^p. 

Additional  direction  cosines  may  be  derived  from  the  Lagrange  Identity 
for  the  expansion  of  the  vector  multiplication, 

(e^  X  I)«  Ce^x  e^)  -  (e^*  ©p)  Ci  •  ®t^ 

As  (e^xl)  -  sin/3pj  and  e^xe^,  •  e^,  equation  (B17)  becomes 
sin)3p  (*cos  21*  12')  *  0  -  cos  (x',r). 


and  so 


cos  (x',r)  •  0.9323  cos)3p. 


The  remaining  direction  cosines  are  derived  from  the  relations 

cos(x',t)cos(x',r)  ♦  cos(y(,t)co3(y',r)  ♦  cos(z,t)cos(z,r)  -  0, 
cos(x',r)cos(y',r)  ♦  cos(x',t)cos(y',t)  ^  cos(x' ,z)cos(y' ,z)  -  0, 
cos(y',r)cos(z',r)  ♦  cos(y',t)cos(z',t)  ♦  cos(y',z)cos{z',z)  *  0, 

and  there  is  obtained 
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cosCa’i^y  «  -  cos  21*12*  -  -  0.9323  sin)9p 
co$(x',z)  -  sin  21*12*  sinP^  -  0,36l6  sin/Sp 
co8(*',a)  ■  ■»  cos/3p  sin  21*12*  ■  -  0,3616  coSySp 

The  above  relationships  may  be  summarized  in  a  table  of  direction  cosines 


r 

1 

t 

z 

X* 

0.9323  3in/3p 

cos  /3p 

0.3616  sinfi 

y' 

0*3616 

0 

-0.9323 

z* 

-0.9323  sin  /9 

P 

8in/?p 

-0.3616  cos  j9p 

For  the  first  pitot-probe  that  was  made,  the  angle  )Sp  -  39*  U6*,  and 
for  the  second,  modified  pitot-probe  there  is  a  choice  of  two  settings 
corresponding  to  angles  for  ^p  of  18*  55'  and  UO*  15', 

The  velocity  components  in  the  (r,t,z)  system  are  calculated  directly 
from  the  equations 


where 


Wr  ■  w^,cos(x*,r)  ♦  Wy,cos(y*,r)  ♦  w^,cos(z',r), 
w^  •  w^,cos(x*,t)  ♦  Wy,,cos(y*,t)  ♦  w^,cos(z' ,t) , 
w^  •  w^,cosCx*,z)  ♦  Wy,cos(y*,z)  ♦  Wg,cos(z*,z), 


(B18) 


w^,  ■  w  cos  0, 

w  ,  ■  w  sin  0  cos§  (BI9) 

y 

w^,  ■  w  sin  0  sin 6* 


To  Illustrate  the  use  of  the  above  analysis,  a  con^>lete  numerical  calcu¬ 


lation  for  one  of  the  data  points  will  be  carried  through. 
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U.  Numerleal  Example 

To  illustrate  the  method  of  computation,  the  flow  data  for  a 
specific  measur'^ment  with  the  pitot-probe  will  be  worked  out  in  detail. 
Consider  the  22CX)  rpm  test-run  having  a  discharge  of  28. li  cfm,  with  a 
pitot  setting  of  r  -  6.720  inch  and  Q  -  The  manometer  read¬ 

ings  are  tabulated  below. 


Manometer  Reading 
(inch  Water-Oage) 

Equivalent  Feet  of 
Air  at  62*F 

h  - 

0 

-  0.28 

-  19.2 

1.90 

130.0 

V'j  ■ 

3.19 

1 

218.3 

■ 

1 

CM 

SS 

0.15 

10.5 

From  Fig.  B5',  for  the  ratios  (h^-h^)/Ch^-h^)  •  0.60  and  (h2-h^)/(h^-h2) 

■  0.05,  the  following  values  are  obtained:  0  •  10.3*>  86.0*  and 

(kg-k^)  ■  1,01.  The  relative  velocity  is  computed  from  equation  (B13), 

w  ■  •  ll8,0  fps. 

The  components  with  respect  to  the  pitot-probe  are  derived  from  equation 
(B19), 

w^,  -  116.1  fps, 

Wy.,  •  1.5  fps, 

w^,  -  21.1  fps. 

The  probe  was  sat  at  the  angle •  U0*l5',  and  from  the  table  of  d^ection 
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eoslnts,  and  aquation  (818),  tho  velocity  conponente  with  respect  to  the 
liqpeller  are 


Wj,  -  55. U  fpe, 

•102.1  fps, 
w.  •  19.9  fpa, 

m 

"fw  w  ■  58.9  fps. 
m  I  r  s 

The  flow  directions  are 


5  •  arc  tan  "  70,2*, 

/9  •  arc  tan  •  29.9*, 

and  the  tangential  component  of  the  absolute  velocity,  and  the  whirl>moment, 


are 


Ct  •  rW  -  w^ 

-  Co, 560) ( 230 .u) ^102.1  •  26.8  fps, 


and 


re. 


15.0  ftVsec. 


The  static  pressure  head,  h  ,  is  determined  from  equation  (515). 

Prom  Fig.  b6,  for  6’  •  86.0*  and  0  •  10,3* >  the  pressure  recovery  factor  is 

k  ■  0.96U.  The  manometer  reading  must  be  corrected  for  the  centrifugal 
0 

effects  upon  the  air  column  in  the  tubing  leading  from  the  pitot  probe  to 
the  rotary  seals.  The  absolute  pressure  at  the  center  hole  of  the  probe  is 


(Po) 


(P. 


)  oxp 


shaft 


gas 


where  R... 
gM 


is  the  gas  constant,  and  the  shaft  radius  is  1.250  inch. 


lOlf 
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67  substitution  into  equation  (Bl^),  and  conversion  to  equivalent  feet  of 
eir«  the  static  pressure  head  is  h  ■  22.2  ft. 

9 

The  quantities  H  and  1  are  obtained  directly  from  the  definitions 

c2 

c  2.W  2 

-  ♦  hg  -  87.2  ft, 

t  .  H  -  i**  m  •  ^20.1  ft. 
g 


I 


I 
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DETAILS  0?  THE  BLADE-TO-BUDE  REUXilTION  SOLUTION 

1.  Method  of  Solution 

The  solution  of  differential  equation  (U3)»  together  with  the  boundary 
conditions  shown  in  Fig.  is  accomplished  by  numerical  procedures  es* 
sentially  developed  by  Southwell  ^  .  The  method  consists  in  estimating 
the  values  of  the  stream  function  over  a  system  of  mesh  points  of  a  grid 
covering  the  flow  region.  The  error  (’’residual"  introduced  into  the 
finite-difference  equation  {hh)  by  these  estimated  values,  is  expressed 
by  the  equation. 

'f^l  *  '1^2  *  '^o  ^  ^^^^0^^  h'  '^3^  "  »  Rg  (Cl) 

The  estimated  values  of  are  ad;}usted  ("relaxed") until  the  residuals  are 
eliminated,  A  grid  spacing,  h  •  0,1,  was  used  Initially  in  the  computation. 
This  grid  spacing  provides  satisfactory  accuracy  in  the  upstream  region 
whore  the  flow  conditions  are  quite  uniform;  however,  in  the  downstream 
regions,  and  at  the  blade  inlet  tip,  the  gradient  of  the  stream  function 
is  relatively  large,  and  in  these  regions  it  is  necessary  to  use  finer 
grid  spacing  for  improved  accuracy.  The  grid  spacing  is  shown  in  Fig,  16. 
At  the  inlet  tip,  and  in  the  downstream  regions,  the  grid  spacing  is  0,0$; 
and  in  the  immediate  neighborhood  of  the  blade  exit  tip,  where  the  Kutta 
condition  must  be  satisfied,  the  grid  spacing  is  0,025» 

2.  Irregular  Grid  Points 

At  the  blade  surface  the  nodes  of  the  relaxation  grid  will  not  coincide 
with  the  blade  contour.  For  these  points,  irregular  star-operators  must 
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be  used,  The  most  general  case,  in  which  all  oi"  the  arms  are  of  non¬ 
standard  length,  is  shown  in  Fig.  Cl.  Considering  the  origin  of  the 
C  C ,  17)  axes  at  the  central  point,  ^  can  be  represented  by  the  polynomial 
approximation, 


*  \  (I)  ♦  ^2 


(C2) 


The  b-coefflclents  are  determined  by  solution  of  the  five  simultaneous 
equations  obtained  from  substitution  of  values  at  the  five  nodal  points. 
This  leads  to  the  expressions  for  the  derivatives. 


y2  2  ^  ^  2  ^ 

u  ^3  4,  *1  0/ _ 

■  a^Xa^Tap’  ^1  ”  a^a^+a”]  ^3  a^a^ 


2 

“  •l‘3 

^0* 

(C3) 

*  ^ 

♦0. 

(CU) 

(a^-a^) 

a<,a. 

% 

(C5) 

For  the  special  case  of  a^  ■  a^  •  1,  the  values  of  the  coefficients  of  ^ 

0 

are  tabulated  by  Shaw°. 

The  blade  inlet  tip' does  not  coincide  with  a  nodal  point  of  the  grid, 
so  the  irregular  star  shown  in  Fig.  C2  is  used.  Applying  the  notation  of 
this  figure  to  equation  (C2),  the  expressions  for  the  derivatives  at  the 
Inlet  tip  become 


h*  7*+  -  ^  V  <-3)  »  ^  [  V  A  \  .(2.V  '*'4 

'*  'll  ■  5*  V  *3’  *  *0’  *  V 
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*1 


T» 


a. 


etc  • 


h  »  standard  arm  length 


0 


Fie  G1  NOTATION  FOR  IREEGUUR 
STAR  OPERATOR 


Tig.  C2 


notation  FDR  STAR 
OPERATOR  APPLIED  AT 
blade  inlet  TIP 


Fig.  C3 


irregular  IRANSITION 
POINT  BETWEEN  FINE  AND 
COARSE  ®IDS. 
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At  transition  points  between  the  fine  and  coarse  grid  spacing,  irreg* 
ular  star  operators  must  be  used  corresponding  to  equations  (G3),(CU),  and 
(C5) .  In  addition,  another  type  of  transition  point  occurs  as  shown  in 
Fig.  C3.  Using  the  polynomial  approximation,  equation  (B2),  the  deriva^ 

tiTSs  are 


4-3-  p  -  ^  4-1) 


(C8) 


(G9) 


h  V  *3>  ‘  h  *ir  *5)  *  ( V‘5»  '''1-  ’''oj 


3.  Kutta-Joukowaky  Condition 


Considering  the  blades  to  be  thin,  the  flow  must  be  tangent  to  the 
blade  surface  at  the  exit.  For  the  preliminary  relaxation  solution,  no  at-- 
tenpt  was  made  to  satisfy  this  condition.  The  downstream  >Aiirl*moraent 


value  was  chosen  on  the  basis  th-t  the  flow  leaves  the  impeller  uniformly, 
at  the  specified  blade  angle.  After  completion  of  the  relaxation  solution, 
it  was  found  that  the  /9 -angle  of  the  streamline,  at  the  blade  exit  tip, 
was  about  10  degrees  less  than  the  blade  angle.  To  correct  this  condition, 
a  decrease  in  the  downstream  whirl-moment  was  necessary.  As  a  first  approx¬ 
imation  of  the  change  in  whirl-moment  required,  consider  a  unit  change  in 
whirl-moment,  and  let  represent  the  corresponding  changes  in  the 
stream  function,  ^^eferring  to  equation  (Cl),  the  values  of  h,  f,  and  g 
are  small,  and  the  terms  involving  these  factors  approach  zero  to  a  higher 
order,  as  the  grid  spacing  is  reduced,  than  the  first  four  terms  of  the 
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equation,  the  flnlte-dlffez*enoe  equation  for  Ai//  then  beeomeB>  approx* 
tnatel7> 

*  ^^2  *  ^'^3  ♦  ®o* 

On  the  blade  contour,  Is  fixed,  ao  that  la  zero  at  the  blade.  Equa¬ 
tion  (Cll) ,  repreaentlng  the  lAplace  equation  in  A  ,  can  be  aolved  by 
the  relaxation  method,  ka  thla  equation  la  linear,  the  aolutlon  for  any 
desired  change  of  idilrl-nonent  can  be  obtained  by  multiplication  by  the  solu¬ 
tion  for  the  unit  change  In  whirl-moment.  After  the  whirl-moment  Is  ad¬ 
justed  to  satisfy  the  Kutta  condition,  a  final  relaxation  must  be  done  to 
eliminate  the  small  residuals  of  eqiiatlon  (Cl). 

U.  Determination  of  the  Velocity  Diatributlon 

Lines  of  constant  ^-values  are  relative  streamlines,  and  the  numerical  f 

value  of  ^  Indicates  the  rate  of  flow  passing  between  the  streamline  and 

* 

the  leading  blade  face,  as  a  percentage  of  the  total  flow  passing  through 
the  channel  between  adjacent  blades.  The  spacing  of  the  streamlines  (see 
Fig.  16)  Is  Indicative  of  the  magnitude  of  the  relative  velocity,  close 
spacing  Indicating  higher  velocities.  The  relative  velocity  components 
along  the  periphery  at  the  blade  outlet  are  determined  by  numerical  differ¬ 
entiation  of  the  stream-function  using  a  seven-point  differentiation  formula 

7 

given  by  Soutluell  for  equally  spaced  grid  points,  and  a  five-point  dlf- 

9 

ferentlatlon  formula  given  by  Wu  for  the  unequally  spaced  grid  points. 
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